NAVAL POSTGRADUATE SCHOOL 
Monterey, California 




THESIS 



THE USE OF ELECTRICAL TRANSMISSION LINE 
THEORY TO PREDICT THE PERFORMANCE 
OF SPACECRAFT RADIATORS 

by 

Steven M. Smith 
March 1992 

Thesis Advisor: Allan D. Kraus 



Approved for public release; distribution is unlimited. 



1254703 



Unclassified 

SECURITY CLASSIFICATION op s page 



REPORT DOCUMENTATION PAGE 



Form Approved 
OMB No 0704-0188 



la REPORT SECURITY CLASSIFICATION 

UNCLASSIFIED 



lb restrictive MARKINGS 



2a SECURITY CLASSIFICATION AUTHORITY 



3 DISTRIBUTION AVAILABILITY OF REPORT 



2b DECLASSIFICATION / DOWNGRADING SCHEDULE 



4 PERFORMING ORGANIZATION REPORT NUMBER(S) 



5 MONITORING ORGANIZATION REPORT NuMBER(S) 



6a NAME OF PERFORMING ORGANIZATION 

Naval Postgraduate School 



6b OFFICE SYMBOL 
(If applicable) 

39 



7a NAME OF MONITORING ORGAN ZATiON 

Naval Postgraduate School 



6c. ADDRESS (City State , and ZIP Code) 

Monterey, CA 93943-5000 



7b ADDRESS (City, State, and ZIPCode) 

Monterey, CA 93943-5000 



8a NAME OF FUNDING SPONSORING 
ORGANIZATION 



8b OFFICE SYMBOL 
(If applicable) 



9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 



8c ADDRESS (City State and ZIP Code) 



10 SOURCE OF FUNDING NUMBERS 



PROGRAM 


PROJECT 


HTa s K 


WORK UNIT 


ELEMENT NO 


NO 


NO 


ACCESSION NO 



11 TITLE (Include Security Classification) 

THE USE OF ELECTRICAL TRANSMISSION LINE THEORY TO PREDICT THE PERFORMANCE OF SPACECRAFT 
RADIATORS (UNCLASSIFIED) 



12 PERSONAL AUTHOR(S) 

Smith, Steven Mark 



13a TYPE OF REPORT 


13b TIME COVERED 


14 DATE OF REPORT (Year, Month Day) 


15 PACE COUNT 


Master's Thesis 


FROM TO 


March 1992 


137 



16 SUPPLEMENTARY NOTATION 

The views expressed in this thesis are those of the author and do not reflect the 
official policy or position of the Department of the Defense or the U.S. Government 



COSAT 1 CODES 


FIELD 


GROUP 


SUB GROUP 















18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 

Thermal Transmission Matrix, Cascade Algorithm, Extended 
Surface, Cooling Fin, Radiative Fin, Longitudinal Fin, 
Rectangular Fin, Trapezoidal Fin, Triangular Fin, Optimum 



19 ABSTRACT (Continue on reverse if necessary and identify by block number) 

The cascade algorithm that is used for extended surface analysis depends on a new parameterization 
called the thermal transmission matrix to represent a single fin. This thermal transmission matrix, which 
is intended to replace the more familiar fin efficiency as a design and analysis parameterization, is a linear 
transformation that maps conditions of heat flow and temperature at the fin tip to heat flow and temperature 
conditions at the fin base. The cascade algorithm was derived by resorting to an analogy between a fin and 
the electrical transmission line. The cascade algorithm permits a fin to be subdivided into many subfins each 
having a thermal transmission matrix and then the individual transmission matrices for each of the subfins 
can be used, via a simple matrix product to form an overall equivalent thermal transmission matrix for the 
entire fin. This thesis develops a thermal transmission matrix for the radiating rectangular, trapezoidal, and 
triangular fins both for the free space and nun-free space environments. Test cases have been run and their 
solutions exactly match those contained in the literature. The thesis concludes with optimization studies for 



20 DISTRIBUTION /AVAILABILITY OF ABSTRACT 

DJjNCLASSIFIED/UNLIMITED □ SAME AS RPT □ OTIC USERS 



21 ABSTRACT SECURITY CLASSIFICATION 

Unclassified 



22a NAME OF RESPONSIBLE INDIVIDUAL 

Allan D. Kraus 



22b TELEPHONE (Include Area Code ) 

(408) 646-2730 



22 c OFFICE SYMBOL 

EC/Ks 



DO Form 1473, JUN 86 



Previous editions are obsolete 

S/N 0102-LF-0 14-6603 



SECURITY CLASSIFICATION OF THIS PAGE 

Unclassified 



Approved for public release; distribution is unlimited. 



The Use of Electrical Transmission Line Theory to Predict 
The Performance of Spacecraft Radiators 

by 

Steven M. Smith 

/' 

Lieutenant, United States Navy 
B.S.E.E., University of Kansas, 1986 

Submitted in partial fulfillment 
of the requirements for the degree of 

MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 
(SPACE SYSTEMS ENGINEERING) 

from the 

NAVAL POSTGRADUATE SCHOOL 
March 1992 



iii 



I MOD 

c.i 



ABSTRACT 

The cascade algorithm that is used for extended surface analysis depends on a 
new parameterization called the thermal transmission matrix to represent a single fin. 
This thermal transmission matrix, which is intended to replace the more familiar fin 
efficiency as a design and analysis parameterization, is a linear transformation that 
maps conditions of heat flow and temperature at the fin tip to heat flow and 
temperature conditions at the fin base. The cascade algorithm was derived by 
resorting to an analogy between a fin and the electrical transmission line. The 
cascade algorithm permits a fin to be subdivided into many subfins each having a 
thermal transmission matrix and then the individual transmission matrices for each of 
the subfins can be used, via a simple matrix product to form an overall equivalent 
thermal transmission matrix for the entire fin. This thesis develops a thermal 
transmission matrix for the radiating rectangular, trapezoidal, and triangular fins both 
for the free space and non-free space environments. Test cases have been run and 
their solutions exactly match those contained in the literature. The thesis concludes 
with optimization studies for each profile considered where it is observed that simple 
algebraic equations can be employed to describe the optimum geometry. 
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I. INTRODUCTION 



Since the advent of the space age, extended surfaces have been extensively 
studied as a means for heat rejection from all sorts of vehicles. Because weight 
optimization is of the utmost importance in spacecraft design, it is desirable to have 
an accurate technique that will rapidly optimize the geometry of the extended surface 
used for rejecting the heat absorbed on or generated within the vehicle. To this point 
in time, the optimization problem involved the time consuming solutions of second 
order non-linear differential equations, by numerical methods or by cumbersome and 
possibly inaccurate graphical methods. A more advantageous method for determining 
the minimum mass, maximum radiant heat rejection extended surface geometry has 
been needed. 

The purpose of this study is to employ the new parameterization called the 
thermal transmission matrix to represent a single fin to evaluate the performance of 
several longitudinal fins. Each fin studied is subdivided into 100 subfins and the 
thermal transmission matrix is calculated for each subfin. Then, by a simple matrix 
multiplication of the individual subfin thermal transmission matrices, the overall 
transmission matrix is computed. Not only will this method produce the base and tip 
conditions of the overall fin entity, but the temperature and heat flow at any point 
between the fin base and fin tip can be determined. This iterative method considers 
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the spacecraft radiant heat rejection problem in both the free and non-free space 
environments. Optimization must be considered and simple algebraic equations that 
describe the fin geometry are formulated for the longitudinal fins of rectangular, 
trapezoidal, and triangular profiles. 

This new approach to fin synthesis and analysis permits an exact solution to any 
degree of accuracy that the user wishes to specify and produces the result in a 
computational efficient manner. 
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H. THEORETICAL PRESENTATION 



A. GENERALIZED DIFFERENTIAL EQUATION 

There are three distinct modes by which heat flows from one point to another. 

A rather simplistic definition of each can be stated. Conduction is the mode by which 
heat flows through a material by molecular motion. If the heat flows from a 
confining surface to a moving fluid, the mode is known as convection. If no 
intervening medium is present and the thermal energy is transported from one point to 
another in the form of electromagnetic waves, the mode is known as radiation or 
radiant heat transfer. 

The flow of heat by conduction is proportional to the temperature gradient and 
cross sectional area perpendicular to the direction of heat flow. 



q cd « - A(x) 



d]W 

dx 



( 2 . 1 ) 



Insertion of a proportionality constant yields 



q cd - - k (2.2) 

where the minus sign assures a positive heat flow in the presence of a temperature 
gradient which must be negative. This serves to define the thermal conductivity of 
the material written as k(T) because it is a material property almost always a function 
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k(T) = 



(2.3) 



Qcd 



- Mx) 



dTQ 

dx 



of temperature. Thus, the heat flow at x normal to a surface of area A(x) at a 
temperature T(x) is written as 



Q C d = ~ k(T) A(x) 



d]\x) 

dx 



(2.4) 



In convection the heat flow from a confining surface at temperature T(x) to a 
fluid whose bulk temperature is T s will be proportional to the temperature difference 
and the amount of surface area. 



« S(x) [ T\x) - T ] 



Insertion of a proportionality constant permits the representation 



Qcr = K six) i m - t s ] 



and this serves to define the convective heat transfer coefficient h cv as 



(2.5) 



( 2 . 6 ) 



^ _ ^£V 

^ " S(x) [ Tlx) - T s ] 



(2.7) 
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and thus Newton’s Law of Cooling is written as 



q cv = hJx,T) S(x) [ 1\x) - T s ] (2.8) 

where the convective heat transfer coefficient h cv is written as a function of position 
and temperature. 

Radiation requires no intervening medium and, for a perfect emitter and 
absorber in full view of one another, the familiar Stefan-Boltzmann law shows that the 
heat How is proportional to the surface area and the difference of the fourth powers of 
the absolute temperature. 

q r = o S(x) [ T, 4 - T 2 4 ] (2.9) 

The proportionality constant used here is the Stefan-Boltzmann constant, 
a = 5.66961 x 10‘ 8 watt / nr K 4 . 

If the surfaces are not perfect emitters and not perfect absorbers and they are 
not in full view of one another, the departure from perfection is via the emissivity 
factor F, and the shape or arrangement factor F A . 

q r = O F a F e Six ) [ r, 4 - r 2 4 ] (2.10) 

The problem of a finned surface radiating to outer space simplifies the foregoing 
equation (2.10). The temperature of the fin becomes T(x) and, because outer space is 
at about 2 K, T 2 can be assumed to be equal to 0. [Compare T(x) 4 = 200 4 to T s 4 = 
2 4 ]. Moreover, F A is equal to 1, the value of a small body in a large enclosure (small 
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Figure 2.1 Energy Balance 

fin in a large enclosure, outer space). Finally the emissivity factor for this case 
(small body in a large enclosure) is F { = e(x,T). The emissivity e is a property of the 
radiating fin surface and temperature. Thus equation (2.10) is rewritten as 

q r = a e(x,7) S(x) T(x) 4 (2-H) 

The difference between the rates of heat entering and leaving the element in 
Figure 2.1 by conduction (dq cdb - dq cda ) must, in the steady state, be equal to the 
rates of heat dissipated from the face, L dx by convection (dq cvl and dq cv2 ) and 
radiation (dq rl and dq r2 ) plus any rates of heat that arrive from external sources (dq sl 
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and dq^). 

The presence of two terms for dq cv , dq r , and dq, permit unequal values on the 
two faces. Thus with 



^cd,a d^cd.b 



dx 



k(x) 6 (x) L 

dx 



dx 



( 2 . 12 ) 



the two convective dissipation terms are 



^cvi + <*?cv 2 = KM* 7 ) Ldx [ T ~ T si] 



+ h„j&T) Ldx[T -T 52 ] 



(2.13) 



and the two radiation dissipations terms with F A = 1 and F t = e(x,T) are 



dq rl + dq r2 = o e(x,T) L dx [ 7 4 -T s ) ] 



+ o e(x,T) L dx 



4 T ^ 



T -T< 



S2 



The energy balance with Tsi — T s2 — 0 is 



(2.14) 



d_ 

dx 



m 5(x) l 

dx 



dx = [ h^ixj) + h^ix, T) ] L dx T 



+ o e(x,T) L dx T 4 + q u + q i2 (2.15) 
In space, if either face of the fin is subjected to an external heat input of the 



form 



q. = a L E ' dx 



(2.16) 
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where E consists of external heat inputs from solar, E s , and/or terrestrial, Ep, sources 
then equation (2.15), the steady state thermal energy balance for a differential element 
dx, is rewritten as 



d_ 

dx 



k(x) b(x) L 



dT 

dx 



dx = 



[ Ki C*.J) 



+ ]L T dx 



+ a e(x,T) LT*dx + 2aLEdx ( 2 - 17 ) 



B. LIMITING ASSUMPTIONS 

Further assumptions must be applied to the generalized differential equation 
(2.17) before a mathematical solution can be attempted. The assumptions employed 
by Murray 1 and Gardner 2 are used for a starting point for further analysis. The 
Murray-Gardner Assumptions listed in Table 1 are essential in simplifying the 
formulation of the generalized differential equation (2.17) so that mathematical 
analysis is possible. Consider an arbitrary symmetrical longitudinal fin shown in 
Figure 2.2 of length (b-a), depth L, and width 5(x) dissipating radiant heat to a 
constant temperature T s surrounding environment. 

The profile area A p is defined as 

A p = J d(x) dx (2.18) 

x=b-a 
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The initial conditions at the fin tip are 



T a = T( x=a ) (2.19) 

and 

q a = q{ x=a ) (2.20) 

Applying the tenth Murray-Gardner assumption (Table 1), that no heat leaves from 
the tip of the fin, equation (2.20) leads to 

dTt - xsa > . 0 ( 2 . 21 ) 

dx 

Two more boundary conditions are taken at the fin base 

7* x=b ) = T b (2.22) 

and 

q{ x=b ) = q b (2.23) 

Applying the fifth and ninth Murray-Gardner assumptions (Table 1), allows one to 
ignore the radiant heat flow from the edges. The surface area S(x) perpendicular to 
the radiant heat flow is given by 

Six) = 2 L y ( 2 - 24 ) 

The face area S(x) through which the heat is dissipated depends on the arc length y, a 
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function of the fin length b. 



The cross sectional area A(x) perpendicular to the heat flow by conduction is 

A(x) = L 6(x) (2-25) 




Because the second Murray-Gardner assumption (Table 1) states that the material is 
homogeneous, both the emissivity and thermal conductivity are assumed to be 
constant and uniform 
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TABLE 1 



THE MURRAY-GARDNER ASSUMPTIONS 4 



1. The heat flow into the fin and the temperature at any point on the fin 
remains constant with time. 

2. The fin material is homogeneous, its thermal conductivity is the same in 
all directions and remains constant. 

3. The heat transfer coefficient between the fin and the surrounding medium 
is uniform and constant over the entire surface of the fin. 

4. The temperature of the medium surrounding the fin is uniform. 

5. The fin width is so small compared with its height that temperature 
gradients across the fin width may be neglected. 

6. The temperature at the base of the fin is uniform. 

7. There are no heat sources within the fin itself. 

8. Heat transfer to or from the fin is proportional to the temperature excess 
between the fin and the surrounding medium. 

9. There is no contact resistance between fins in the configuration or between 
the fin at the base of the configuration and the prime surface. 

10. The heat transferred through the outermost edge of the fin (the fin tip) 
is negligible compared to that through the lateral surfaces (faces) of the fin. 
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e(x,T) = e 



(2.26) 



k(T) = k (2.21) 

The contribution from convective dissipation from the faces of the fin is ignored 
because the arbitrary symmetrical longitudinal fin is operating in a space environment. 

C. GENERALIZED DIFFERENTIAL EQUATION FOR A FIN 

Combining the generalized differential equation (2.17) and equations (2.24) 
through (2.27), and assuming that there are no external heat inputs (qj, = q^ = 0), 
the steady state thermal energy balance for the differential element shown in Figure 
2.2 is 

—Ik L 6(*)^^ -2 Lae T(x) 4 = 0 (2.28) 

dx\ dx ) dx 

The temperature excess 6(x) is defined as the temperature difference between the 

temperature at point x and the constant temperature of the surrounding 

environment T s . 

0(x) = T(x) - T s (2.29) 

and from this, it is seen that 

dl\x) = dQ(x) (2-30) 

Applying the temperature excess equations (2.29) and (2.30), the equation for the 
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steady state thermal energy balance for the differential element is rewritten as 



- Uh. aw 4i = o (2.3i) 

k dx 

where h r is defined as the radiative heat transfer coefficient. Because T s = 0, 

h r = o e 6(x) 3 (2.32) 

A fundamental concept of calculus is that the length of arc d 7 in a cartesian 
coordinate framework is 



A 6(I) 

dx dx 



dy = [ dx 2 + dy z ] 1/2 
The derivative of the arc length is 



it X , lim^l = ± 1 + 

dx ax-o Ax \ \dx) 



(2.33) 



(2.34) 



The + or - sign is to be taken according as the arc length y increases or decreases as 
x increases. 

Combining the equation for the steady state thermal energy balance differential 
element (2.31) and the arc length equation (2.34), the generalized differential equation 
for a longitudinal fin of arbitrary symmetrical profile is 



d_ 

dx 



b(x) 



dQ(x) 

dx 



2 h 



- 0(x) [l + &(x) 2 ] m = 0 



(2.35) 



Equation (2.35) is a second order, linear, homogeneous differential equation with the 
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boundary conditions 



and 



rf9( x = a > = 0 (2.36) 

dx 



0( x=b ) = e b (2.37) 



D. LONGITUDINAL FIN OF RECTANGULAR PROFILE 

The longitudinal fin of rectangular profile is illustrated in Figure 2.3. It should 
be noted that the origin of the x axis is located at the tip of the fin and not at the base 
of the fin. Because the width is constant over the entire profile, the width function 
5(x) is simply 

5(x) = 6 b = 6 a (2.38) 

The profile area A p is written as 

x=i 

A p = f 5(x) dx = 5 b b ( 2 -39) 

x=0 



Inserting the width function (2.38) into the generalized differential equation for a fin 
equation (2.35) yields 



* dx 2 



IK 

k 



0(x) = 0 



(2.40) 
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or 



where 



- WJ 2 0(JC ) = o (2.41) 

dx 2 



m 




(2.42) 



The radiative heat transfer coefficient h r is assumed to be a constant for a very small 
subfin where T b ~ T a . 

The general solution of the generalized differential equation (2.41) is 



0(x) = C x e mx + C 2 e'" 1 * (2-43) 

where Q and C 2 are arbitrary integration constants that are evaluated by applying the 
initial condition of equation (2.22). 



d(b) = Q b = C x e mb + C 2 e' mb (2.44) 

The heat flow, q b , into the base of the fin is found by inserting the general 
solution of the generalized differential equation (2.43) and the cross sectional area 
perpendicular to the heat flow by conduction (2.25) into Fourier’s Law of 
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Figure 2.3 Longitudinal Fin of Rectangular Profile 



Conduction, equation (2.4). 

q(x) = L b b k | -^[Cj e mx + C 2 e' m *]j (2.45) 

Simplifying and evaluating at x = b with the initial condition of equation (2.23), 
equation (2.45) becomes 

q b = L b b k m ( C, e mb - C 2 e' m b ) (2.46) 
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Solving for C, and C 2 yields 



C, 



, -m b 



0*. + 



-m b 



2 Y, 






(2.47) 



and 



where 



«m b 



c 2 - 



0. - 



,m b 



2 Y. 



<h 



(2.48) 



Y b = L ( 2 h r k 6 b ) 1/2 (2.49) 

is defined as the characteristic thermal admittance Y b of the fin. Likewise the 
characteristic thermal impedance, 7^, of the fin is the reciprocal of Y b 

z ‘ = i . = i ( 2 K k > ' m (2 - 50) 

An examination of equations (2.43) through (2.48) reveals that only two of the 
four variables are independent and that the specification of any two variables allows 
for the determination of the remaining two. Combining equations (2.43) through 
(2.48) yields the two port inverse thermal transmission matrix for a longitudinal fin of 
rectangular profile 



0- 




cosh mb 


-Z b sinh mb 




% 


4a 




-Y b sinh mb 


cosh mb 




q b. 



(2.51) 
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where the inverse thermal transmission matrix is defined as 



T" 1 = 



cosh mb 
-Y b sinli mb 



-Z b sinh mb 
cosh mb 



(2.52) 



The thermal transmission matrix is defined as the inverse of the inverse thermal 



transmission matrix. 



T 



= T" 1 = 



A 

C 



B 

D 



(2.53) 



where A, B, C, and D are the elements of the thermal transmission matrix. The 
thermal transmission matrix representation is 



e; 




cosh mb 


Z b sinh mb 






q b. 




Y b sinh mb 


cosh mb 




% 



(2.54) 



E. LONGITUDINAL FIN OF TRAPEZOIDAL PROFILE 

The longitudinal fin of trapezoidal profile is shown in Figure 2.4. Note that the 
origin of the axis is not at the tip of fin but located to the right of the tip 
where x = 0. 

The width function <5(x) for this profile is 

6(jc) = 2 x tan 4> (2.55) 
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where </> is defined as 



4> = arctan 



*> 

2 (b - a) 



(2.56) 



The profile area A is 



x-b 

A p = f b(x) dx = ~ a) {b b + bj 



(2.57) 



Inserting the width function (2.55) into the generalized differential equation for a fin 
[Equation (2.35)] yields 



or 



where 



2 x tan(4>) 



J 2 6(x) 

dx 1 



+ 2 tan(4>) 



dm 

dx 



— ©c*) - 0 (2.58) 

k cos(4>) 




dm 

dx 



- m 2 0(x) = 0 



(2.59) 



m 



1/2 



k sin(4)) 



(2.60) 



Equation (2.59) represents a second order, linear, homogenous differential equation. 
The general solution is expressed in terms of modified Bessel functions. The general 
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Figure 2.4 Longitudinal Fin of Trapezoidal Profile 5 



solution to equation (2.59) is 

0(x) = Cj I 0 { 2 m x 1/2 ) + C 2 K 0 { 2 m x 1/2 ) (2.61) 

The following transformation is made 

u = 2 m x 1/2 (2.62) 

such that at x = b 

u b = 2mb^ (2-63) 
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and at x = a 



u a = 2 m a 1/2 (2.64) 

After applying the boundary conditions of equations (2.36) and (2.37) and 
making the transformations of equations (2.63) and (2.64), the general solutions 
become 

= C, W ♦ C 2 W (2.65) 

and 

2 k b.L m 2 

<7* = [C, /,(«*) - C 2 *,(«*)] (2-66) 

u b 



Simultaneous solution of the general solution equations (2.65) and (2.66) yield 



^1 U b 



W 0, 



2 k 6. L m 2 

< O V 



W <i b 



(2.67) 



and 



^2 



W 9, - 



U L 



2 k h.L m 2 

. O i 



K («*) <i b 



(2.68) 



Substituting and rearranging the general solution equations (2.65) and (2.66) and the 
equations for C, and C 2 (2.67) and (2.68), the elements of the thermal transmission 
matrix for the longitudinal fin of trapezoidal profile are written as: 
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A - u a [ /,(«.) KJjiJ * />,) K x {u a )] 



(2.69) 



B = 



2 k b b L m 2 



[WW - W W] ( 2 - 70 > 



C = 2 k bL m 2 

a. 




[ I x {uJ K x {u a ) - I x (U a ) K x ( Ub )] (2.71) 



*>-«*[ W W + W 



(2.72) 



F. LONGITUDINAL FIN OF TRIANGULAR PROFILE 

The longitudinal fin of triangular profile shown in Figure 2.5 is a special case of 
the longitudinal fin of trapezoidal profile. There are many similarities between the 
longitudinal fin of triangular and trapezoidal profiles. The differential 
equation for the longitudinal fin of triangular profile is the same as the differential 
equation for the longitudinal fin of trapezoidal profile equation (2.59). 




dx 2 dx 



(2.73) 
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Figure 2.5 Longitudinal Fin of Triangular Profile 6 



where 



m 



k sin(4>) 



1/2 



(2.74) 



and 



4 > 



arctan 



A. 

2 b 



(2.75) 
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The profile area A p is 



x=b 

A p = f 6(x) dx = | 6 b 

x=0 Z 



(2.76) 



However the evaluation of the arbitrary constants from the boundary conditions 
differ between the triangular and the trapezoidal profiles. As the value of x 
approaches zero the temperature excess at the tip 0 a must be equal to some finite 
value, therefore 

I 0 ( u= 0 ) = 1 (2.77) 

In order to physically obtain a finite temperature at x a = u a = 0, C 2 in equation 
(2.61) must be equal to zero because K 2 (0) is unbounded. The unbounded condition 
allows cancellation of the K c (u) terms from the general solutions for temperature 
excess and heat flow. 

After cancellation and rearranging the general solutions for the longitudinal fin 
of trapezoidal profile equations (2.65) and (2.66), the general solution for the 
longitudinal fin of triangular profile becomes 



Vb = 



2 k h b L m 2 



W 

W 



(2.78) 



As a result of the zero tip width, the edge parameters do not depend on the base 
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parameters. The thermal input admittance element for the longitudinal fin of 
triangular profile is 



Y = Si = 
' e t 



2 k b b L m 2 7,(up 
u b 7 o ( u i) 



(2.79) 
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m. CASCADE ANALYSIS ALGORITHM 



A. GENERALIZED TWO PORT TRANSMISSION LINE MODEL 

The radiating fin is a distributed heat flow configuration which is analogous to 
the transmission line which is a distributed electrical network. The general two port 
transmission line model is shown in Figure 3.1. The model is described by two input 
terminals. At point "a" there is an applied voltage V, and a current I 4 flows into the 



- -A. 




-A - 




Transmission 








V b '•» 




Line 











Figure 3.1 Two Port Transmission Line Model 
model. There is also a pair of terminals at point "b" where a voltage V b occurs and a 
current I b flows out of the terminals. The two port model of the transmission line can 
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consist of any arrangement of passive elements within the two port but with no active 
current or voltage sources within the two port. 

It can be shown from network theory that the two port transmission line model 
in Figure 3.1 can be represented by the following equations: 

V a = A V b + B I b (3.1) 



I a = C v b + D I b 

where the determinant of the ABCD matrix defined by 







A 


B 




v b 


A 




C 


D 




A 



is equal to unity. 



A D - B C = 1 



(3.2) 



(3.3) 



(3.4) 



The values A, B, C, and D are constants and are defined as the elements of the 
transmission matrix. 

The electrical transmission line analogy can be directly applied to an equivalent 
thermal transmission line model where the heat flow is analogous to current and the 
temperature is analogous to voltage. Equations (3.1) and (3.2) are rewritten in the 
equivalent thermal transmission line model as: 

T a =■ A T b + B q b (3.5) 
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q a = C T b + D q b 



( 3 . 6 ) 



B. CASCADE ALGORITHM 

The thermal transmission parameters are very useful in describing two port 
networks connected in a cascaded arrangement. It is an exercise in elementary two 
port theory. Examine the two equations for the two networks shown in Figure 3.2. 



and 



CD 

Jr 






K 


i 

CD 

Q> 


q K 




c i 


P 

i 


i 

•Q 

Q) 







a 2 


B 2 


■ 0 a; 






C 2 


p 
1 





( 3 . 7 ) 



( 3 . 8 ) 



Continuity for the configuration in Figure 3.2 dictates that 



and 



which can be represented as 



** = 



( 3 . 10 ) 




( 3 . 11 ) 



The overall thermal transmission parameter matrix for the two networks 
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cascaded in Figure 3.2 is by definition 



A 


B 




X 


X 


to 


X 


C 


D 




X 


X 


,C 2 


X 




Figure 3.2 Cascaded Network 

and the thermal transmission matrix for a cascaded network of n networks is 

A B 
C D 

Equation (3.13) is an expression of what is referred to as the cascade algorithm. 

C. COMPUTER ANALYSIS 

The cascade algorithm expression (3.13) can be applied to a fin divided into n 
subfins each having the same length. Converting the temperature excess notation 









^2 ^2 




X B n 




X X 
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back to the temperature notation (T s = 0), the cascade algorithm expression (3.13) is 
written as 







CQ 

H 


a 2 b 2 




*n B n 


T a 


Qb 




Ci D, 


C 2 D 2 




cj c 


Qa. 



The base temperature T b and the base heat flow q b are the conditions at the base and 
the tip temperature T a and tip heat flow q a are the conditions at the edge. 

A problem arises because the tip temperature T a and the base heat flow q b are 
not known or not specified. The cascade algorithm expression (3.14) can not be 
solved directly because it is essentially a linear transformation that maps tip conditions 
to base conditions. However, a solution can be calculated indirectly. If one makes 
an initial estimation for the tip temperature T a , it can be used to calculate the radiant 
heat flow q a out of the tip of the fin. 

Q a = o e 5 a L T a 4 (3.15) 

Thus T a and q a are the starting values for the computations using the cascade 
algorithm expression (3.14). At the conclusion of the process, a comparison 
involving all n subfins is then made between the known or specified base temperature 
T b and the calculated base temperature. If the comparison indicates an error, the tip 
temperature T a is changed and the iterative process is repeated and continues until the 
error between the presumed known base temperature T b and the calculated base 
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temperature are within an acceptable range. Figure 3.4 provides a complete flow 
chart of the cascade algorithm. 

Further amplification concerning the radiative heat transfer coefficient needs to 
be made. Recalling equation (2.32) where 

h r = o e T(x ) 3 ( 3 . 16 ) 

it is not clear as to what temperature, T(x), should be applied to the computation of h, 
and hence to the cascade algorithm expression (3.14). The left and right faces of the 
differential element in Figure 3.3 have an average temperature. 

T + T. 

T = —2 * ( 3 . 17 ) 

av 2 




Figure 3.3 Differential Element 
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and as long as T a and T b are not vastly different, equation (3.16) can be written using 



T 



h z = o e Tj (3.18) 

The environmental parameter V is introduced for the case in which there are 
portions of the surrounding environment at some temperature other than T s = 0. An 
example of this is in the interception of thermal radiation from a nearby solar panel at 
some temperature not near absolute zero. The environmental parameter is defined as 



V = 



*2 

2 a e T b * 



(3.19) 



where the parameter K 2 allows for the application of external heat inputs from solar, 
E s , or terrestrial, E T , sources or both. 



K 2 = E a 



(3.20) 



The environmental parameter is equal to zero for the free space case and to some 
number greater than zero but less than one for the non free space case. Substituting 
equation (3.19) into (3.18), the radiant heat transfer coefficient becomes 



h z - a e 



*2 

2 T av 



(3.21) 



External heat inputs (q u and q^) from solar input E s and terrestrial input Ej- are now 
applied to the cascade algorithm expression (3.14) through the application of the 
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environmental parameter V. 



D. PROGRAM DESCRIPTION 

The computer code consists of three analysis programs for each of the profiles 
considered, and a MAIN Program that is used for general housekeeping. The 
PERFORMANCE Analysis Program evaluates an entered set of parameters returning 
the characteristics (volume, profile area A p , profile number PN, base heat flow q b , 
edge temperature T a , and the fin efficiency rj) of the fin along with the thermal 
transmission parameter matrix. The PROFILE AREA Analysis Program also 
evaluates an entered set of parameters but it returns a data set of 100 fins having a 
constant profile area A p with the base fin width 5 and length b which are used as the 
input variables. The OPTIMIZE Analysis Program is similar to the 
PERFORMANCE Analysis Program with the exception that it returns the optimum 
fin geometry which results in the minimum mass and maximum base heat flow fin. 

The MAIN Program provides menu driven access to the various analysis 
programs. The user has the ability to set the computational main program error 
coefficient. The user also has the ability to set the save data switch which controls 
data flow into an ASCII formatted disk file. Commas are used for the data delimiter 
which allows file importation into MATLAB, Word Perfect, Draw Perfect, and Math 
CAD for additional analysis or for plotting. An additional feature provided is the unit 
conversion module. Table 2 list the available units that are provided to the user by 
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the MAIN Program. Figure 3.4 provides a flow chart of the cascade algorithm 
expression (3. 14) which is common to all of the analysis programs. 

The PERFORMANCE Analysis Program requires the following input conditions 
to be entered: 

• base temperature 

• density 

• emissivity 

• environmental parameter 

• fin base width 

• fin edge width 1 

• fin depth 

• fin length 

• thermal conductivity 

The program calculates and returns the following values: 

• base heat flow 

• edge temperature 

• fin efficiency 

• profile area 

• profile number 

• thermal transmission parameter matrix 

• total mass 

The fin efficiency r\ and the profile number PN are extensively referred to in 
literature and are provided for use in making comparisons. The fin efficiency and the 
profile number are defined as 



<jb 

2 o e b L T b * 



(3.22) 



'Applies for only the longitudinal fin trapezoidal profile. 
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Figure 3.4 Cascade Algorithm Flow Chart 
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TABLE 2 



UNITS AVAILABLE IN THE MAIN PROGRAM 



Admittance 


• watts / °C, BTU / hr °F 


Area 


• m 2 , cm 2 , ft 2 , in 2 , mm 2 


Heat 


• watts, kw, BTU / sec, BTU / hr 


Impedance 


• °C / watts, °F hr / BTU 


Length and Depth 


• m, cm, ft, in, mm 


Mass 


• kg, grams, lbm, oz 


Temperature 


• K, °C, °F, °R 


Thermal Conductivity 


• watts / m K, BTU / ft hr °F 


Width 


• m, cm, ft, in, mm 



PN 



2 O G b 2 T b * 

1 



(trapezoidal) 



(3.23) 



PN = 






o 6 -b 2 Tt, 3 (rectangular) (3.24) 
k b b 



If the save switch is switched on, the following differential element conditions are 
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written to a disk file: 



• base heat flow 

• base temperature 

• distance from the tip of the fin 

The OPTIMIZE Analysis Program requires the following input conditions to be 
entered: 



• base temperature 

• density 

• emissivity 

• environmental parameter 

• fin edge width 2 

• fin width 

• thermal conductivity 

The program calculates and returns the following optimum values: 

• base heat flow 

• edge temperature 

• fin base width 

• fin efficiency 

• fin length 

• profile area 

• profile number 

• thermal transmission parameter matrix 

• total mass 

The program saves the same quantities that the PERFORMANCE Analysis Program 
saves depending on the condition of the save switch. 

The PROFILE AREA Analysis Program requires the following input conditions 



2 Applies for only the longitudinal fin trapezoidal profile. 
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to be entered: 



• base temperature 

• density 

• emissivity 

• environmental parameter 

• profile area 

• starting fin base width value 

• ending fin base width value 

• taper ratio 3 

• fin depth 

• thermal conductivity 

The program calculates and writes to a disk file the following parameters for 100 fins 
having the same profile area Ap within the specified width range: 

• base heat flow 

• base width 

• base temperature 

• edge temperature 

• fin efficiency 

• thermal transmission parameter matrix 

• total mass 

The analysis programs all use the thermal transmission matrices computed in 
Chapter 2 in the cascade algorithm with the exception of the longitudinal fin of 
triangular profile. The triangular and trapezoidal profile analysis are very similar. 

The analysis program uses the thermal admittance equation (2.79) to compute the base 
conditions for the first subfin having a thickness of The subfin is located at the tip 
of the fin. The base conditions are then applied and used as the tip conditions in the 
thermal transmission matrix for the trapezoidal profile (2.69) through (2.72). The 



3 Applies for only longitudinal fin of trapezoidal profile. 
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cascade algorithm is then used to calculate the thermal transmission matrices for the 
remaining 99 subfins in order to finally determine the overall thermal transmission 
matrix and the base conditions. The trapezoid that is being used in the calculation has 
an almost zero tip width compared to the base width such that 

6 b = 100 6 a (3.25) 

Borland Turbo C++ was used to write the computer code. The computer code 
listing appears in Appendix A. Full compiler optimization was selected in order to 
decrease the run times. Average run times on an IBM 486/487 machine appear in 
Table 3. 



TABLE 3 

AVERAGE PROGRAM RUN TIMES 





PERFORMANCE 

(sec) 


PROFILE AREA 
(sec) 


OPTIMIZE 

(sec) 


Rectangular 


4 


97 


4 


Trapezoidal 


4 


600 


4 


Triangular 


5 


625 


5 
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IV. PROGRAM ANALYSIS RESULTS AND DISCUSSION 



A. PROGRAM ERROR TOLERANCE 

In order to allow for meaningful interpretation, an understanding of the accuracy 
of the results produced from the cascade algorithm is vital. A test to show 
convergence is used to illustrate the approximate relative error. 

Table 4 provides an illustration of how the error coefficient in the MAIN 
Program affects the accuracy of the program results. Table 4 is produced by 

TABLE 4 

PROGRAM ERROR TOLERANCE 



Program 



Error 


Ta 


* % 


q b 


« % 


Coeff % 


(K) 


error 


(watts) 


error 


1.000 


282.18 


0.0106 


926.28 


2.0918 


0.500 


282.24 


0.0106 


931.83 


1.5052 


0.100 


282.19 


0.0071 


946.36 


0.0307 


0.050 


282.22 


0.0035 


946.94 


0.0920 


0.010 


282.22 


0.0035 


946.10 


0.0032 


0.005 


282.21 


0.0000 


946.05 


0.0021 


0.001 


282.21 


0.0000 


946.07 


0.0000 



running the same input parameters many times while decreasing the error coefficient 
in the MAIN Program. The PERFORMANCE Analysis Program for the longitudinal 



52 



fin of rectangular profile was used. The Table 4 results are applicable for any of the 
analysis programs because of the commonality of the cascade algorithm. As the error 
coefficient in the MAIN Program decreases, the output value converges. The 
approximate relative error e t is defined as 



C r s 



True Value - Approximate^ Value 
True Value 



100 



(4.1) 



The exact true value is unknown, however the value produced when the error 
coefficient in the MAIN Program is equal to 0.001 can be assumed to be 
approximately equal to the exact true value. 

Examining the results illustrated in Table 4, one can conclude that when the 
error coefficient in the MAIN Program is less than or equal to 0.1 the approximate 
relative error observed in the output is generally much less than the error coefficient 
in the main program. 



B. GRAPHICAL EXAMPLES 

The purpose of this section is to describe and provide examples on how each of 
the three analysis programs can be used. The computer code listings are provided in 
Appendix A. 

1. Example One: Longitudinal Fin of Rectangular Profile Analysis 

A longitudinal fin of rectangular profile is 0.5 meter long, 1.0 centimeter 
thick, and 1.0 meter deep. The fin is made of aluminum (k = 209.4 watt / m K, 
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p = 2.6 kg / dm 3 , and e = 0.85). If the base temperature is 400 K, what are the tip, 
midpoint, and base conditions (error < 0.1%) both in the free space and non free 
space environments? 

Because the geometry of the fin is already known, the PERFORMANCE 
Analysis Program is used for the evaluation. The save switch in the MAIN Program 
must be placed in the save position in order to later retrieve the midpoint conditions 
from the disk file. The computer output display for the free space environment 
(V=0) appears in Table 5. 

Figures 4. 1 and 4.2 provide a detailed graphical representation of the 
required conditions both in the free space and non free space environments. The base 
conditions are read from where the length is equal to 50 centimeters and the tip 
conditions are read from where the length is equal to 0 centimeter. In the free space 
environment the tip, midpoint, and base conditions for this particular fin geometry are 
listed in Table 6. 
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TABLE 5 



PERFORMANCE ANALYSIS OUTPUT DISPLAY 



THE LONGITUDINAL FIN 


OF 


RECTANGULAR PROFILE 


Density 


= 


2.60 kg / dnT3 


Fin Depth 


= 


1.00 m 


Total Mass 


= 


13.00 kg 


Fin Length 


= 


0.50 m 


Fin Width 


= 


1.00 cm 


Emissivity 


= 


0.85 


Profile Area 


= 


50.00 cm~2 


Profile Number 


= 


0.607 


Base Heat Flow 


= 


732.85 watts 


Fin Efficiency 


= 


0.594 


Base Temperature 


= 


400.07 K 


Edge Temperature 


= 


324.45 K 


Thermal Conductivity 


= 


209.40 watts / m K 


Environmental Parameter 


= 


0.00 


Transmission Parameter Matrix 


A = 1.2287 


B 


i = 0.2583 deg C / watt 


B = 2.2374 watt / deg C 


D 


i = 1.2842 



TABLE 6 

FIN CONDITIONS IN FREE SPACE 



Location 


Temp 


q 




(K) 


(watts) 


Base 


400.07 


732.85 


Midpoint 


341.66 


292.33 


Tip 


324.45 


10.68 
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Figure 4.1 Heat Dissipation From A Rectangular Profile Fin 



56 




(51 2sp) du»x 



Figure 4.2 Temperature Distribution For A Rectangular Profile Fin 
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Equation (3.4) is used below to verify that the determinant of thermal 



transmission matrix listed in Table 5 is indeed equal to unity. 



det 



1.2287 0.2583 

2.2374 1.2842 



1.0000 



(4.2) 



2. Example Two: Longitudinal Fin of Rectangular Profile Area Analysis 
A longitudinal fin of rectangular profile 1.0 meter deep is made of 
aluminum (k = 209.4 watt / m K, p = 2.6 kg / dm 3 , and e = 0.85). It is required 
to remove a minimum of 500 watts of heat from the base of a fin at a temperature of 
400 K. Determine the geometry of the fin using several values for the profile area 

A P . 



The PROFILE AREA Analysis Program is used because the geometry of 
the fin is unknown. Figure 4.3 provides a graphical representation which shows the 
constant profile area A p relationship as a function of width versus heat flow. Any of 
the dimensions listed in Table 7 would meet the minimum requirement to remove the 
500 watts of heat from the base of the fin. 

Examining Figure 4.3, one can see that the maximum heat flow occurs at 
a particular width value. More heat is dissipated as the profile area Ap is increased. 
The value of the width at the point of maximum heat flow also increases as 
the profile area A p increases. The point of maximum heat flow is considered as the 
optimum geometry (maximum heat flow, minimum mass) for the particular profile 
area A p . A summary of the optimum geometry points extrapolated from the data files 
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FIGURE 4.3 Profile Area Analysis 
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TABLE 7 



PROFILE AREA ANALYSIS POINTS 



A P 


s b 


b 


q b 


(cm 2 ) 


(cm) 


(m) 


(watts) 


25 


0.334 


0.748 


507.4 


100 


0.490 


2.041 


633.7 


200 


0.724 


2.762 


771.2 


300 


0.334 


8.982 


523.1 



produced from the PROFILE AREA Analysis Program appear in Table 8. 

TABLE 8 

OPTIMUM GEOMETRY POINTS FOR THE RECTANGULAR PROFILE 



A P 


^tx> 


b 0 


Mass 0 


q bo 


qjm 


(cm 2 ) 


(cm) 


(m) 


(kg) 


(watts) 


(w/kg) 


25 


0.646 


0.387 


6.50 


579.76 


89.19 


100 


1.543 


0.648 


26.00 


924.53 


35.56 


200 


2.440 


0.819 


52.00 


1165.83 


22.42 


300 


3.103 


0.967 


80.00 


1334.40 


16.68 



3. Example Three: Comparison of Various Profiles 



A longitudinal fin 1.0 meter deep having a profile area A p equal to 50 cm 2 
is made of aluminum (k = 209.4 watt / m K, p = 2.6 kg / dm 3 , and 
e = 0.85). The base temperature is held at 400 K. Compare the relative 
performance between the longitudinal fins of rectangular, trapezoidal, and triangular profiles. 
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The same approach is used as in the previous example. A longitudinal fin 
of trapezoidal profile having a taper ratio TR = 0.5 is used. The taper ratio is 
defined as 

TR - — (4.3) 

b b 

Figure 4.4 provides a graphical representation which shows the relationship between 
the three profiles being considered. 

The curve in Figure 4.4 produced by the longitudinal fin of trapezoidal 
profile is midway between the curves produced by the longitudinal fins of rectangular 
and triangular profiles. Of more importance is the fact that at the location of the 
optimum geometry points, the longitudinal fin of triangular profile produces the 
largest heat flow while the longitudinal fin of rectangular profile produces the smallest 
heat flow. The optimum geometry points extrapolated from the disk file generated 
from the PROFILE AREA Analysis Program appear in Table 9. Not only is there an 
optimum fin geometry for a particular profile area A p but there is also an optimum 
profile shape which produces maximum heat flow. 

4. Example Four: Optimum Geometry Comparison 

A longitudinal fin 2.0 meters deep (e = 0.95) is required to remove at 
least 1000 watts from a spacecraft in a free space environment (V=0). The base 
temperature is held at 400 K. Using the list of materials appearing in Table 10 
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Figure 4.4 Comparisons of Various Profiles 
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Width (cm) 



TABLE 9 



OPTIMUM GEOMETRY POINTS FOR TIIE VARIOUS PROFILES 



Profile 

RECT 

TRAP 

TRI 



Taper 

Ratio 

1.0 

0.5 

0.0 



S b 

(cm) 

0.958 

1.270 

1.309 



b 

(m) 

0.522 

0.524 

0.763 



732.47 

772.93 

826.78 



q b /m 

(W/kg) 

56.343 

59.456 

63.598 



Mass q b 
(kg) (watts) 

13.00 
13.00 
13.00 



conduct a heat flow per unit mass comparison using the longitudinal fins of 
rectangular and triangular profiles. 



TABLE 10 

MATERIALS USED IN THE OPTIMUM COMPARISON 8 





Chemical 


P 


k 


Material 


Symbol 


kg/dm 3 


w / m K 


Aluminum 


A1 


2.6 


209.4 


Copper 


Cu 


8.8 


384.0 


Iron 


Fe 


7.2 


58.0 


Magnesium 


Mg 


1.7 


157.0 


Molybdenum 


Mo 


10.2 


145.0 


Silicon 


Si 


2.3 


83.0 


Silver 


Ag 


10.5 


407.0 


Zinc 


Zn 


6.8 


140.0 
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The OPTIMIZE Analysis Program is used for the analysis. The 
parameters were entered individually into the OPTIMIZE Analysis Program and the 
results are listed in Table 1 1 . 



TABLE 11 

OPTIMUM COMPARISON RESULTS 







*bo 




Mass 0 


q b o 


qbo/m 


Material 


cm 


cm 2 


kg 


watts 


w / kg 


Ag 


RECT 


0.40 


21.73 


45.63 


997.86 


21.86 


Ag 


TRI 


0.54 


15.09 


31.69 


995.45 


31.41 


A1 


RECT 


0.79 


42 . 24 


21.96 


998 .96 


45.49 


A1 


TRI 


1.05 


29.33 


15.25 


995.47 


65.27 


CU 


RECT 


0.43 


23 . 03 


40.54 


997 .94 


24 . 61 


Cu 


TRI 


0.57 


15.99 


28.15 


995.45 


35.36 


Fe 


RECT 


2 . 84 


152 .49 


221.11 


1003 . 99 


4.54 


Fe 


TRI 


3.79 


105.89 


153 . 54 


995.75 


6.48 


Mg 


RECT 


1.05 


56.33 


19.60 


999.74 


51.01 


Mg 


TRI 


1.40 


39.12 


13.61 


995.49 


73 . 14 


Mo 


RECT 


1.14 


61.00 


124.43 


1000.00 


8 . 04 


Mo 


TRI 


1.51 


42 .36 


86.41 


995.54 


11.52 


Si 


RECT 


1.99 


106.56 


49.66 


998.41 


20.11 


si 


TRI 


2.65 


73.99 


34.48 


995.59 


28.87 


Zn 


RECT 


1.18 


63 . 17 


85.92 


1000.12 


11.64 


Zn 


TRI 


1.57 


43.87 


59.66 


995.50 


16.68 
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Figure 4.5 shows a bar graph representation of the heat flow per unit mass 
relationship using the various materials listed in Table 10. It is evident from Figure 
4.5 that the optimized longitudinal fin of triangular profile made of magnesium is the 
most efficient material considered of those listed in Table 10, with aluminum being a 
close second. The melting temperature of both materials (*= 900 K) is far above the 
base temperature (400 K) stated in the problem statement. Comparing Tables 10 and 
11 materials having a small value for the density (p < 3.0 kg / dm 3 ) and a large 
value for the thermal conductivity ( k > 150 watt / m K) produce large values for 
heat flow per unit mass. 

C. DERIVATION OF WORKING OPTIMUM GEOMETRY EQUATIONS 

As shown in Tables 8 and 9 optimum geometry points can be determined using 
the PROFILE AREA Analysis Program. The overall success of finding the 
dimensions for the optimum geometry using the PROFILE AREA Analysis Program 
has several limitations. The first limitation is that the starting and ending fin base 
width values must be known. And it can be argued that a second limitation is that 
excessive analysis program run times which might exceed ten minutes. 

Relationships for the optimum geometry for the longitudinal fins of rectangular, 
trapezoidal, and triangular profiles would remove both limitations. Kern and Kraus 9 
and Chung and Nguyen 10 stated that the dimensions for optimum fin geometry can 
be expressed by 



65 




Figure 4.5 Optimum Profile Comparison Using Various Materials 
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where Pi(TR), P 2 (TR), and P 3 (TR) are defined as polynomials in which the taper ratio 
TR is the independent variable. 

In general the polynomial P is expressed in the form of 



P(TR) = Cj TR + C 2 TR " ~ 2 + ••• + C, 77? 0 (4.7) 

where n is the degree of the polynomial. 

The polynomials P,, P 2 , and P 3 are determined from running the PROFILE 
AREA Analysis Program using a slightly modified code. The computer code listing 
appears in Appendix B. A random number generator is used to generate the input 
parameters. The random number generator produces input parameters in the range 
listed in Table 12. Murray-Gardner (Table 1) assumption number five (the fin width 
is small compared with its depth and length) is checked in order to insure that the 
original assumptions are kept valid. 
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TABLE 12 



RANGE OF THE RANDOM INPUT PARAMETERS 



0.01 < e < 0.99 
20 < L ( cm ) < 220 
10 < k ( watt / m K) < 210 
100 < T b ( K ) < 600 
25 < q b ( watts ) < 1025 
0.01 < Taper Ratio TR < 0.99 



The program enters into a loop calculating the heat flow q b using the random 
input parameters. The program remains in this loop until the maximum heat flow q b 
is found. Once found the program writes the random input parameters and the value 
determined to be the maximum heat flow q b to a disk file. Another set of random 
input parameters is generated and the maximum heat flow q b once again is calculated. 
This sequence continues until 100 optimum data sets are written to a disk file. 

Three data sets using the random input parameters are produced, one each for the 
longitudinal fins of rectangular, trapezoidal, and triangular profiles. 

The last step in determining the polynomials Pj, P 2 , and P 3 is to fit the collected 
data stored in the disk file to a "best" polynomial fit having the smallest degree d 
using the least square principle. Equations (4.4) through (4.6) are rewritten as 
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(4.8) 



P(TR) 



L bo e T. A 
0 0 

% 




(4.9) 




(4.10) 



Equations (4.8) through (4.10) are placed into a MATLAB "m" file. Each of the 
collected data files are read into MATLAB and were evaluated using the MATLAB 
polynomial fitting function polyfit with the taper ratio TR as the independent variable 
using varying degrees d in order to find the polynomials P,, P 2 , and P 3 which "best" 
approximate the previously collected data from the disk file. The results appear in 
Table 13. The listing of the MATLAB evaluation m files appear in Appendix C. 

The standard error e s is defined as 



1/2 



E Oi - y .) 2 



(4.11) 



N - 1 



where y ; is the computed value of yj and there are N values. 
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TABLE 13 



COEFFICIENTS OF THE OPTIMUM POLYNOMIALS 





Degree 








Standard 


Profile 


d 


c, 


c 2 


c 3 


Error 


PI RECT 


2 


0.8675 


0.0000 


0.0000 


3.54e-3 


PI TRAP 


2 


-0.0290 


0.8138 


0.0000 


2. 12e-3 


PI TRI 


2 


0.0000 


0.9040 


0.0000 


2.55e-3 


P2 RECT 


2 


1.8648 


0.0000 


0.0000 


1.63e-2 


P2 TRAP 


3 


0.3881 


-1.0324 


2.7522 


1.22e-2 


P2 TRI 


2 


0.0000 


2.4852 


0.0000 


1.40e-2 


P3 RECT 


2 


1.6178 


0.0000 


0.0000 


2.06e-2 


P3 TRAP 


2 


0.5147 


1.1382 


0.0000 


1.24e-2 


P3 TRI 


2 


0.0000 


1.1234 


0.0000 


9.46e-3 



The equation for the optimum profile area A^ for either the rectangular, 
triangular, or trapezoidal profiles is written as 



= ( 1.1382 + 0.5147 TR ) 



( \3 


1 


% 




k o 2 e 2 


L T\ 

\ b ) 





(4.12) 
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V. CONCLUSIONS 



This thesis has described a new technique that uses a parameterization called the 
thermal transmission matrix to represent a single fin to evaluate the performance of 
several fins. The algorithm is based on the simple ideal of subdividing the fin into 
100 subfins and then the thermal transmission matrix is calculated for each subfin. 
Then, by a simple matrix multiplication of the individual subfin thermal transmission 
matrices, the overall thermal transmission matrix is computed. Not only will this 
produce the base and tip conditions of the overall fin entity, but the temperature 
distribution and heat flow at any point between the fin base and the fin tip can be 
easily determined. 

Several conclusions can be made from the Chapter 4 Program Analysis Results 
and Discussion. It was observed that the profile area Ap, increases as the cube of the 
heat flow q b . In order to double the heat flow q b , the profile area A p must be 
increased eight times. 

Material selection is very important. The weight of the fin is proportional to 
the density and the thermal conductivity is inversely proportional to the profile area. 
Materials having a thermal conductivity greater than 150 watt / m K and a density 
less than 3 kg / dm 3 are generally the best suited materials. For the same material, 
environmental conditions, heat flow, and base temperature, the optimum longitudinal 
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fin of triangular profile requires only 69 percent as much material as the optimum 
longitudinal fin of rectangular profile, with the optimum longitudinal fin of 
trapezoidal profile falling in between. 

Simple working relationships to compute the optimum fin dimensions are found. 
Table 16 compares the findings from this thesis with the results from previous work 
by others. The differences between the optimum dimension polynomial coefficients 
listed in Table 16 are very small, which verifies the accuracy of the thermal 
transmission matrix parameterization technique. 

TABLE 16 

OPTIMUM DIMENSION POLYNOMIAL COEFFICIENTS 



RECTANGULAR 



TRIANGULAR 



Pi P2 P3 



.8675 


1.8648 


1.6178 


This Thesis 


.8846 


1.8473 


1.6341 


Chung and Nguyen 1 


.8842 


1.8579 


1.6056 


Liu 12 


.8846 


1.8473 


1.6341 


Bartas and Sellers 13 


.9040 


2.4852 


1.1234 


This Thesis 


.9485 


2.3104 


1.0957 


Chung and Nguyen 


.9546 


2.3027 


1.0991 


Nilson and Curry 14 


.9091 


2.4255 


1.1025 


Kern and Kraus 15 



This new approach to fin synthesis and analysis permits an exact solution to any 
degree of accuracy that the user wishes to specify and produces the result in a 
computational efficient manner. 
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APPENDIX A COMPUTER CODE LISTING 



j ★ ***★******★*************★**★***★***********★**★*★★*★*★★★*★*★★**★★*★★★★★★★★* 



Main Program 

A*************************************************************************** * j 



#include <process.h> 
^include <conio.h> 
#include <stdlib.h> 
^include <stdio.h> 

^define Yes 1 
^define No 0 



char 


*S1 □ 


= 


c 


II 


watt / deg C", » BTU / hr deg F" ); 




char 


*S2[] 


= 


c 


II 


m A 2" 


, " cm A 2", " ft A 2", " in A 


2", " n*n A 2" ) 


char 


*S3 [] 


= 


c 


II 


watts" 


, " kw", " BTU / sec", " BTU / 


hr" }; 


char 


*S4[] 


= 


c 


II 


deg C 


/ watt", " deg F hr / BTU" }; 




char 


*S5 [] 


= 


i 


II 


m", " 


cm", " ft", " in", " mm" }; 




char 


*S6 [] 


= 


t 


II 


kg", " 


grams", " Ibm", " oz" ); 




char 


*S7 [] 


= 


t 


II 


deg K" 


, " deg C", " deg F", " deg R" 


); 


char 


*S8 [] 


= 


t 


II 


watt / 


m deg K", " BTU / ft hr deg F" 


1 ); 


char 


*S9 [] 


= 


t 


II 


No", " 


Yes"); 




char 


*0[] 


= 


i 




"<1> 


watt / deg C <2> BTU / hr deg 


F", 



"<1> m A 2 <2> cm A 2 <3> ft A 2 <4> in A 2 <5> mm A 2", 

"<1> watts <2> kw <3> BTU / sec <4> BTU / hr", 

"<1> deg C / watt <2> deg F hr / BTU", 

"<1> m <2> cm <3> ft <4> in <5> mm", 

"<1> kg <2> grams <3> Ibm <4> oz", 

"<1> deg K <2> deg C <3> deg F <4> deg R", 

"<1> watt / m deg K <2> BTU / ft hr deg F", 

"<1> No <2> Yes" ) ; 

void print_main_menu( ) , print_menu_one( ) , pri nt_menu_two( ); 

void setupO, overlaysO, updatejjni ts( ); 

float error=0.1; 

int T=0, L=0, M=0 , S=0, H=0, d=4, K=0, A=1, 2=0, Y=0; 

char main_key=0; 

int main( K 

FILE *data_in; 

if ( (data_in = f open( "uni ts . ext" , "r" ) ) != NULLK 

f scanf (data_in, "%d %d %d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,&d,&K,&A,&Z,&Y,&error); 
f c lose(data_in) ; ) 

pr i nt_ma i n_menu( ) ; 

do{ 

V if (main_key == '1') setupO; 

if (main_key >= '2' && main_key <= '4') overlaysO; ) 

while (( main_key = getchO ) != '5'); 

windowd ,1,80,25); 

cl rscr( ) ; 

_setcursortype(_NORMALCURSOR ) ; 
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return 0;} 



void over lays( void){ 
char key=0; 

int Quit = No, result; 
pr i nt_menu_one( ) ; 
do 

switch (key){ 



case 'A 1 : cl rscr( ); 
if (main_key 
if (main_key 
if (main_key 
break; 



Quit = 


Yes; 


r= '2 1 ) 


result 


== 1 3 * ) 


result 


== >4') 


result 



= spawn l (P_WAIT, 
= spawn l (P_WAIT , 
= spawnl (P_WAIT, 



"RECT.EXE", NULL); 
"TRAP.EXE", NULL); 
"TRI.EXE", NULL); 



case 1 B 1 : clrscrO; Quit = Yes; 



if 


(main_key == 1 2 ' ) 


result 


if 


(main_key == 1 3 1 ) 


result 


if 


(main key == '4 1 ) 


result 


break; 





spawn l (P_WAIT # "ARECT.EXE", NULL); 
spawn l (P_WAIT, "ATRAP.EXE", NULL); 
spawn l (P_WAIT , "ATRI.EXE", NULL); 



case 'C': clrscrO; 
if (main_key 
if (main_key 
if Cmain_key 



Quit = Yes; 

== '2') result 
r= *3') result 
== '4') result 



spawnl (P_WAIT, "ORECT . EXE" , NULL); 
spawn l (P_WAIT , "0TRAP.EXE", NULL); 
spawnl (P_WAIT, "0TRI.EXE", NULL);) 



while (( key = toupper(getch( ) ) ) != 1 D 1 && Quit == No); 



print_main_menu( ); 
return;) 

void setup( void){ 
char key=0; 



pr i nt_menu_t wo C ) ; 
do 

switch (key){ 



case 1 A 1 : 



case ' B 1 : 



case 1 0: 



case 1 D 1 : 



case • E • : 



case 'F‘: 



case ' G 1 : 



gotoxy(2, 17); clreolO; cprintf ("%s",0[0] ); 
while (key != '1' && key != '2') key = getch(); 

Y = key - 1 1 1 ;print_menu_two( ); break; 

gotoxy( 1 , 17); clreolO; cprintf ("%s",0[1] ); 

while (key != *1' && key != 1 2 ' && key != '3'&& key != '4 •&& key != '5') 
key = getch( ); 

A = key - ' 1 ' ;print_menu_two( ); break; 
gotoxy(1 , 17); clreolO; cprintf ("%s",0[2] ); 

while (key != 'V && key != '2'&& key != *3 •&& key != '4') key = getchO; 
H = key - 1 1 ' ;print_menu_two( ); break; 

gotoxy(2, 17); clreolO; cprintf ("%s",0 [3] ); 
while (key != '1* && key != '2') key = getchO; 

Z = key - ' 1 1 ;pr int_menu_two( ) ; break; 

gotoxy(2, 17); clreolO; cprintf ("%s",0[4] ); 

while (key != '1* && key != '2'&S key != *3'&& key != • 4 ' && key != '5') 
key = getchO; 

L = key - 1 1 * ;print_menu_two( ) ; break; 
gotoxy(3, 17); clreolO; cprintf ("%s",0[5) ); 

while (key != '1* && key != '2'&S key != '3‘&& key 1= '4') key = getchO; 
M = key - ' 1 1 ;print_menu_two( ); break; 

gotoxyO , 17); clreolO; cprintf ("%s",0 [6] ); 

while (key != '1' &S key != '2'&& key != '3 •&& key != '4') key = getchO; 



case 1 H 1 : 



case ' 1 1 : 



case 1 J ' : 



case '<•: 



T = key - ' 1 ' ;print_menu_two( ); break; 

gotoxyC 1 , 17); clreolO; cprint*f ("%s",0 [7] ); 
while (key != '1* && key != '2') key = getchO; 

K = key - 1 1 1 ; print_menu_two( ) ; break; 

gotoxyC 2, 17); clreolO; cprintf ("%s",0 [4] ); 

while (key != 1 1* && key != ‘Z'&S, key != 1 3 1 && key != '4'&& key != '5') 
key = getchO; 

d = key - ' 1 ' ;print_menu_two( ); break; 

gotoxyC 12, 17); clreolO; cprintf ("%s",0[8] ); 
while (key != '1' && key != *2') key = getchO; 

S = key - ' 1 ' ;print_menu_two( ); break; 

gotoxy(22, 12); clreolO; scant ("%f",&error); 
print_menu_two( );} 



while (( key = toupper( getchO ) ) 1= 1 L 1 ) ; 



update_uni ts(); 
print_main_menu( ); 
return;} 



void updatejjni ts( void)( 

FILE *uni t ; 

unit = fopen("uni ts.ext", "w"); 

fprintf (unit,"%d %d %d %d %d %d %d %d %d %d %f",T,L,M,S,H,d,K,A,Z,Y,error); 

fclose(uni t); 

return;} 

void print jnain_menu( vo id){ 

_setcursortype(_NOCURSOR); 
windowO , 1,80,25); 
cl rscr( ); 

gotoxy(27,3); cprintf ("EXTENDED SURFACE ANALYSIS"); 
wi ndow( 20,6,80,25); 



cprintf ("<1> 
cprintf ("<2> 
cpr i ntf ("<3> 
cpri ntf ("<4> 
cprintf ("<5> 
return;} 



Setup / Ini tial izat ion\r\n"); 

Fin of Rectangular Prof i le\r\n" ) ; 
Fin of Trapezoidal Prof i le\r\n") ; 
Fin of Triangular Prof i le\r\n") ; 
Return to DOS"); 



void pr int_menu_one(void){ 

cl rscr( ); 
gotoxyC 1,1); 



cprintf (" \r\n" ) ; 

cprintfC" <A> Predict the Perf ormance\r\n" ) ; 

cprintf (" <B> Analyze the Profile Area\r\n"); 

cprintfC" < C> Optimize the Dimensions\r\n" ) ; 
cprintfC" <D> Return to Main Menu\r\n"); 

cprintfC" "); 

return;} 



void pr int_menu_two(void){ 



cl rscr(); 
gotoxyC 1,1); 
cprintfC" — 
cprintf ("<A> 
cpri ntf ("<B> 
cprintf ("<C> 
cprintf ("<D> 
cprintf ("<E> 



Admittance 

Area 

Heat 

Impedance 

Depth 



: %s\r\n", SI [Y] ); 
: %s\r\n", S2[A]); 
: %s\r\n", S3[H]); 
: %s\r\n", S4[Z]); 
: %s\r\n" , S5[L]); 



\r\n"); 



75 



cprintf ( M <F> 
cprintf ("<G> 
cprintf ( M <H> 
cpr i ntf ( "<I> 
cprintf ("<J> 
cprintf ("<K> 
cprintf ("<L> 
cprintf (" — 
return;} 



Hass : %s\r\n", S6[M]); 

Temperature : %s\r\n", S7[T]); 

Thermal K : %s\r\n" f S8[K]); 

Width : %s\r\n», S5[d]); 

Save Data : %s\r\n", S9[S]); 

Percent Error: %f\r\n" f error); 

Return to Main Menu\r\n M ); 



") 
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Jk A************************************************************************** 



PERFORMANCE Analysis Program 
Longitudinal Fin of Rectangular Profile 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk k j 



U include <stdio.h> 
^include <stdlib.h> 
^include <math.h> 
^include <conio.h> 
#include “units.h" 



#def ine 
#def ine 
#def ine 
#def ine 



Yes 1 

No 0 

N 100 

sigma 5. 66961 e- 8 



char f i le_name[20] ; 

double A_()7 B_(), C_(), D_(), m(), hr(), Yb(), A, B, C, D; 

double length, width, depth; 

double inc, overall_diff = 100, Tmax = 0, Tbvjnin = 100000, data [N+1] [6] ; 

double qb=0,Al, B1, Cl, D1, Ta, Tb, Ttip, Tave, temp, qa, segment_di f f ; 

float error, V, epsilon, k, Tbase, density, mass, K2, eta, PN; 

int T, L, M, $, H, d, K, i, X, 2, Y; 

int overflow, possible_overf low = No, number_of_tries=0; 

void inputO, computet ), resultsO; 



Jk kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



MAIN PROGRAM 



*★★****★★★*****★★★***★*******★★★★********★★********★★****★*★★★★******★★**★ * j 



int main(K 

FILE *data_in, *out; 

data_in = f open("uni ts.ext*', "r"); 

fscanf (data_in, u %d %d %d %d %d %d %d %d %d %d %f M ,&T ,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error) 

fclose(data_in); 

error = error / 100; 

i nput ( ) ; 

K2 = 2 * epsilon * sigma * V * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if ( V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/ 2); 

else 

Ttip = .0001 + Tave; 
inc = .5 * Ttip; i = 1; 

while ( fabs(overal l_dif f ) > (Tbase * error) \\ i <= N K 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * width * epsilon * pow(Ta, 4) * depth; 

while (i <= N && overflow == NoK 

segment_diff = 10; 

while (fabs(segment_di f f ) > ((Tb * error) / N) && overflow == No){ 

A = A_( ); B = B_( ) ; C = C_(); D = D_(); 
temp =A*Ta+B* qa; 
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qb =C*Ta+D* qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segmented i f f ; 

if (Tb > Tbase + 50 j | Tb < OK 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbvjnin ) Tbvjnin = Ttip; 
inc = (Tbvjnin - Tmax) / 2; 

Ttip = fabs( Tbvjnin - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == NoK 

data [ i ] [OD = qb; data [ i D [ID = Tb; data[i][2] = A; 
data [ i D [3D = B; data [i] [4] = C; data [i] [5] = D; 
Ta = Tb; qa = qb; ++i;» 



if (number_of_t ries > 29 K 

clrscrO; cprintf ("Fatal Error: Unable to Compute"); 

getch(); exit(0);> 

if (overflow == No)( 

if (Tmax < Ttip) 

if ( Tb < Tbase ) Tmax = Ttip; 

overal l jdiff = (Tb - Tbase) / 2; 

if ( possible_overf low == YesK 
inc = inc / 2; 

if (overal IjJif f >0) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overal l_di ff; 

number_of _t r i es++ ; 

if (Ttip < Tave) Ttip = Tave; 

gotoxy(27 # 24) ; 

cpri ntf ("N = %d Ta = %4.1f Tb = %4.1f " , number_of_tr i es, w_T(T, T tip), w_T ( T , Tb) ) ;}} 

for (i = N - 1; i >= 1; i - - K 

A1 = A * data [i ] [2] + B * data [i ] [43 ; B1 = A * data[i][3] + B * data[i][5]; 

Cl = C * data [i ] [2] + D * data[i] [43; D1 = C * data[i] [3] + D * data[i] [53 ; 

A = A 1 ; B = B 1 ; C=C1; D=D1;> 

if (S == IK 

out = f open( f i le_name, "w"); 
for (i = 1; i <= N ; i++) 

fprintf(out, "%e, %e, %e\n", w_L(d, i * length / N) ( w_T(T,data [i] [1] ), w_H(H, 

data [i 3 [0] ) ); 

fclose(out);} 

Ta = Ttip; 

mass = width * length * depth * density * 1000; 

PN = length * sqrt( sigma * epsilon * pow(Tbase,3) / ( k * width)); 
eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase # 4 ) ); 
resul ts( ); 
exi non- 
return; } 

j-k ★***★★***★★★*******★★★********★★*★★★*★****★*★★★**★★*★★★★★*★**★*★****★★★ 

OUTPUT PROCEDURE 

a ************************************************************************* * j 



void resul ts( vo id)C 



cl rscr( ); 

cprintf ("\a THE LONGITUDINAL FIN OF RECTANGULAR PROFILE"); 

cpr i nt f ( "\r\n ============= =====================================================\r\n" ) ; 
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cpr i nt f C *' Density = %4.2f kg / dm A 3\r\n",density); 

cprintfC" Fin Depth = "); p_LCL, depth); 

cprintfC" Total Mass = "); p_M(M, mass); 

cprintf (" Fin Length = "); p_LCL, length); 

cprintfC" Fin Width = "); pJ_Cd, width); 

cprintfC" Emissivity = %4.2f\r\n", epsilon); 

cprintfC" Profile Area = "); p_ACX, width* length); 

cprintfC" Profile Number = %4.3f\r\n",PN); 

cprintfC" Base Heat Flow = "); p_HCH, qb); 

cprintfC" Fin Efficiency = %4 .3f \r\n", eta) ; 

cprintfC" Base Temperature = "); p_TCT,Tb); 

cprintfC" Edge Temperature = "); p_TCT,Ta); 

cprintfC" Thermal Conductivity = "); p_KCK, k); 

cprintfC" Environmental Parameter = %4.2f\r\n", V); 

cprintfC" ==================================================================\r\n" ) • 

cprintfC" Transmission Parameter Matrix \r\n\r\n"); 

cprintfC" A = %4.4f B = ", A); p_ZCZ,B); 

cprintfC" C = "); p_Y(Y,C); 

cprintfC" D = %4.4f", D); 

getchc ); 

return;) 



I it ************************************************************* ********** 

INPUT PROCEDURE 

************************************************************************** * j 



void input CvoidK 

clrscrC ) ;_set cursor typeC_NORMALCURSOR); 



cprintf C"\r\n THE LONGITUDINAL FIN OF RECTANGULAR PROF I LE\n\r" ) ; 

cprintf C"\r\n ==============================^===================================\r\n" ) ; 

cprintfC" Fin Depth "); depth = i_LCL); 

cprintfC" Fin Length "); length = i _L CL); 

cprintfC" Fin Width "); width = i _L C d ) ; 

cprintfC" Emissivity = "); scanfC"%f", Sepsilon); 

cprintfC" Density C kg / dm A 3 ) = "); scanfC'^f", &density); 

cprintfC" Base Temperature "); Tbase = i_T C T ) ; 

cprintfC" Thermal Conductivity "); k = i_KCK); 

cprintfC" Environmental Parameter CO to .8) = "); scanfC"%f", &V); 

if CS == 1) C 

cprintfC" Output Filename is = "); 



scanfC"%s", file_name);) 



cprintf C"\ r\n==================================================================") ; 

if C epsilon > 1.0 \ \ epsilon <= OK 

cpri ntf C "\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 

getchc ); exitCO);) 

if C V < 0 | | V > 0.8K 

cprintf C"\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 

getchc ) ; exitCO);) 

_setcursortype(_NOCURSOR ); 
return;) 

y* *********************************************************************** 

GLOBAL FUNCTIONS 

************************************************************************** * j 



double hrCvoidK 

float Tav = CTb + Ta) / 2; 

returnCsigma * epsilon * powCTav, 3) - C K2 / C2 * Tav)));) 
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double m(void)( 

return( sqrt( ( 2.0 * hr()) / ( k * width) ) );) 



double Yb(voidK 

return( sqrt( 2.0 * hr() * k 

double A_(voidK 

return( cosh( m( ) * length / 

double B_(voidK 

return ( sinh( m( ) * length / 

double C_(voidK 

return ( sinh( m( ) * length / 

double D_(voidK 

return( cosh( m( ) * length / 



* width) * depth );) 
N ) );> 

N ) / Yb() );> 

N ) * Yb() );> 

N ) );) 
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j-k *************************************************************************** 



PERFORMANCE Analysis Program 
Longitudinal Fin of Trapezoidal Profile 

******************** A******************************************************* it j 



^include <stdio.h> 
^include <stdlib.h> 
^include <math.h> 
tfinclude <conio.h> 
^include "bessel .h" 
# include "units.h" 



#def ine 
#def ine 
#def i ne 
#def ine 



Yes 1 

No 0 

N 100 

sigma 5.66961e-8 



char f i le_name [20] ; 

double A_(), B_(), C_(), D_(), m() # hr(), uCdouble x), A, B, C, D; 

double length, depth, edge_width, base_width, fin_taper; 

double bl , db, a, b, ua, ub, Ta, Tb, a_width, b_width; 

double qb=0, inc, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 100000, data[N+1][6] 

double A1, Bl, Cl, D1, taper_ratio, Ttip, temp, qa, segment_di f f ; 

float V, epsilon, k, error, Tbase, density, mass, K2, Ap, eta, PN; 

int T, L, M, S, H, d, K, i, X, Z, Y; 

int overflow, possi ble_overf low = No, number_of_tri es=0; 

void input (), resultsO; 



jit ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★it********'*************** 



MAIN PROGRAM 



•k-k-kitit'kitititititit'k'k'kitit-kit-kit-k-k-kit-k-k'k-k'k'k'kit'k'k-k-kifkifkifkifkickick'k'k'kititickit'kickicititititickitifkifkitit it j 



int main()( 

FILE *data_in, *out; 

data_in = fopen("units.ext","r"); 

fscanf (data_in, "%d %d %d %d %d %d %d %d %d %d Xf'^&T.&L.&M.&S.&H.&d.&K.&X.&Z.&Y.&error); 

fclose(data_in); 

error = error / 100; 

input (); 

taper_ratio = edge_width / base_width; 
bl = length / ( 1 ~ taper_ratio ); 
fin_taper = atan2( base_width, 2*b1 ); 
db = length / N; 

K2 = 2 * epsilon * sigma * V * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if ( V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/ 2); 

else 

Ttip = T ave + 0.0001 ; 
inc = .5 * Ttip; i=1; 

while (fabs(overal l _d iff ) > (Tbase * error) | | i <= N )( 



overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * edge_width * epsilon * pow(Ta, 4) * length; 

a = bl - length; 
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while (i <= N && overflow == No)( 



segment_diff = 10; 

a_width = edge_width * a / ( bl - length); 

b_width = edge_width * ( a + db ) / ( bl - length); 

while (fabs(segment_diff ) > ( Tb * error) / N && overflow == NoK 

ua = u(a) ; ub = u(a + db); 

A = A_( ); B = B_( ); C = C_(); D = D_(); 
temp = A * Ta + B * qa; 

qb = C * Ta + D * qa; 

segment_diff = ( temp - Tb) / 2; 

Tb = temp + segmentjdi f f ; 

if (Tb > Tbase +50 \\ Tb < 0)( 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbvjnin ) Tbv_min = Ttip; 
inc = (Tbvjnin - Tmax) / 2; 

Ttip = fabs( Tbvjnin - inc); 
if (Ttip < TaveT Ttip = Tave;}} 



if (overflow == No)( 

data [i] [0] = qb; data Ci ] [13 = Tb; data[i] [2] = A; 
data [ i 3 [33 = B; data [ i ] [43 = C; data[i] [5] = D; 
Ta = Tb; qa = qb; a = a + db; ++i;» 



if ( number j>f_tries > 29 X 

clrscrO; cprintf ("Fatal Error: Unable to Compute"); 

getchO; exit(0);} 

if (overflow == NoX 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overf low == YesX 
inc = inc / 2; 

if (overal l_di f f > 0 ) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overal l _d iff; 
number_of _t r i es++ ; 
gotoxy(27,24) ; 

cprintf ("N = %d Ta = %4.1f Tb = %4 - 1 f " # number_of_tries,w_T(T,Tt ip),w_T(T,Tb)) 
if (Ttip < Tave) Ttip = Tave;}} 

for (i = N - 1; i >= 1; i--)< 

A1 = A * data [i ] [2] + B * datali] [4]; B1 = A * data[i] [3] + B * datali] [5] ; 

Cl = C * datali] [2] + D * datali] [4]; D1 = C * datali] [3] + D * datali] [5] ; 

A = A1 ; B = B1 ; C=C1; D=D1;} 

if (S == IX 

out = f open( f i le_name, "w"); 
for (i = 1; i <= N; i++) 

fprintf (out , "%e, %e # %e\n", w L(d, i * length / N), w_T(T , data Ci] [1] ) , w_H(H, 

datali] 10] ) ); 

f close (out);} 



Ta = Ttip; 

Ap = .5 * length * ( edge_width + base_width ); 
mass = Ap * length * density * 1000; 

PN = 2 * sigma * epsilon * pow( length, 2) * pow(Tbase,3) / (k* base_width); 
eta = qb / ( 2 * sigma * epsilon * length * length * pow(Tbase,4)); 
resul ts( ); 
ex i t ( 0 ) ; 
return;} 
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OUTPUT PROCEDURE 



************************************************************************** t, j 



void results(void){ 



clrscr( ) ; 
cprintfC 
cprintf ( M \r\n 
cprintf (" 
cprintf (" 
cpr int f (" 
cprintf (" 
cprintf (" 
cpr int f (" 
cprintf (" 
cprintf (" 
cpr int f (" 
cprintf (" 
cprintf (" 
cprintf (" 
cpr int f (" 
cpr intf (" 
cprintf (" 
cprintf (" 
cprintf (" 
cprintf (" 
cprintf (" 
cprintf (" 
getch( ) ; 
return; > 



THE LONGITUDINAL FIN OF TRAPEZOIDAL PROFILE"); 



Density = %4. 
Fin Depth = "); 
Total Hass = "); 
Fin Length = "); 
Emissivity = %4. 
Prof i le Area = ") ; 
Profile Number = %4. 
Base Heat Flow = ") 
Fin Efficiency = %4. 
F in Edge Width = " ); 
Fin Base Width = "); 
Base Temperature = "); 
Edge Temperature = "); 
Thermal Conductivity = ") 
Environmental Parameter = %4. 



2f kg / dnr3\r\n",densi ty); 
p_L ( L , depth); 
p_M ( M , mass); 
p_L ( L , length); 

2f\r\n", epsilon); 
p_A(X # Ap) ; 

3f\r\n",PN); 

P_H( H , qb); 

3f\r\n",eta); 
p_L(d, edge_width); 
p_L(d, base_width); 
p__T (T # Tb) ; 
p_T (T # Ta); 
p_K(K, k); 

2f \r\n", V); 



Transmission Parameter Matrix \r\n\r\n"); 
A = %4.4f B = ", A); p__Z(Z,B); 

C = "); p_Y ( Y , C ) ; 

D = %4.4f", D); 



=\r\n" ) ; 



/* 



*★★****■* 



********^ 



INPUT PROCEDURE 



******** 



************* 



**************** 



*********** 



void input(void)( 

clrscr( ); _se t cursor type(_NORMAL CURSOR); 



cprintf ("\r\n THE LONGITUDINAL FIN OF TRAPEZOIDAL PROF I LE\r\n" ) ; 

cpr intf ("\r\n ===========r====r=r===rr======r=========================r==r======\r\n" ) 

cpr int f(" Fin Depth "); depth = i_L(L); 

cprintfC" Fin Length "); length = i_LCL); 

cprintfC Fin Edge Width "); edge_width = i_L(d); 

cprintfC Fin Base Width "); base_width = i_L(d); 

cprintfC Emissivity = "); scanf("%f", &epsilon); 

cprintfC Density ( kg / dm A 3 ) = "); scanf("%f", &density); 

cprintfC Base Temperature "); Tbase = i_T(T); 

cprintfC Thermal Conductivity "); k = iJ((K); 

cprintfC Environmental Parameter (0 to .8) = "); scanf("%f", &V); 

if CS == 1) C 

cprintfC Output Filename is = "); 



scanf("%s", file_name);> 



cprintf C\r\n =================================r====================r==s==s=====») ; 

if ( edge_width >= base_width)( 

cprintf C\r\n\r\nFATAL ERROR: Fin Edge Width >= Fin Base Width"); 

getch(); exit(0);) 

if ( epsilon > 1.0 [j epsilon <= 0){ 

cprintf C\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 

getchC); exit(O);} 
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if ( V < 0 ] | V > 0.8X 

cprintf ("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8") 

getchC); exit(0);} 

__setcursortype(_NOCURSOR ) ; 
return;} 

j* *********************************************************************** 

GLOBAL FUNCTIONS 

************************************************************************** * j 



double hr(voidX 

float Tav = (Tb + Ta) / 2; 

return(sigma * epsilon * powCTav, 3) - (K2 / (2 * Tav)));} 
double m(voidH 

returnC hr() / (k * sin(f in_taper)) );} 

double u( double x)( 

returnC 2 * sqrtCx * m( ) ) );} 

double A_(void)( 

returnC ua * ( 1 1 (ua)*K0(ub) + I0(ub)*K1 (ua) ) );} 

double B_(void)C 

return C ub*ub*( I0Cub)*K0(ua) - l0(ua)*K0(ub) ) / ( 2*k*b_width*depth*m( ) ));} 
double C_(void){ 

return ( 2*k*a_width*depth*mC )*ub*( 1 1 (ub)*K1 (ua) - 1 1 (ua)*K1 (ub) )/(ua) );} 
double D_(void)C 

returnC ub*( l0Cua)*K1 Cub)+I 1 (ub)*K0Cua) ) );} 
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Jk kkkkkkkkkkkkkkkkkkkkkkkkkk*kkkkkkk*kkkkk*kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



PERFORMANCE Analysis Program 
Longitudinal Fin of Triangular Profile 

**************************************************************************** * j 



#include <stdio.h> 

# include <stdl ib.h> 
tfinclude <math.h> 
#include <conio.h> 
#include "bessel.h" 
# include "units.h" 



#def ine 
#def i ne 
#def ine 
#def ine 



Yes 1 

No 0 

N 99 

sigma 5. 66961 e- 8 



char f i le_name[203 ; 

double A_( ) , B_(), C_(), D_(), m( ) , hr(), u(double x), A, B, C, D; 

double length, depth, edge_width, base_width, fin_taper; 

double bl, db, a, ua, ub, Ta, Tb, a_width, b_width; 

double qb=0, inc, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 100000, data[N+1][6] 

double A1 , Bl, Cl, D1, T t f p , temp, qa, segment_di f f ; 

float V, epsilon, k, error, Tbase, density, mass, K2, Ap, PN, eta; 

int T, L, M, S, H, d, K, i, X, Z, Y; 

int overflow, possible_overf low = No, number_of_tr i es=0; 

void input(), resultsO; 



jk kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



MAIN PROGRAM 



************************************************************************** * j 



int main(K 

FILE *data_in, *out; 

data_in = fopenO'units.ext'V'r"); 

f scanf (data_i n, "%d %d % d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,&d,&<,&X,&Z,&Y,&error); 
fclose( data-in) ; 
error = error / 100; 

i nput ( ) ; 

db = length / ( N + 1 ); 

bl = length; 

length = length - db; 

fin_taper = atan2( base_width, 2*b1 ); 

K2 = 2 * epsilon * sigma * V * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if ( V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = Tave + 0.0001; 
inc = .5 * Ttip; i=1; 

while ( f abs(overal l_di f f ) > ( Tbase * error) | [ i <= N X 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 
segment_diff = 10; 

while ( fabs(segment_di f f ) > ( Tb * error) / N && overflow == NoX 
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ub = u(db); 

temp = Ta / 10(ub); 

C = 2 * k * edge_width * depth * pow(m(),2) * I1(ub) / (ub * 10(ub)); 
qb = C * temp; 

segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_di f f ; 

if (Tb > Tbase +50 \\ Tb < OK 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbv_min ) Tbvjnin = Ttip; 
l'nc = (Tbvjnin - Tmax) / 2; 

Ttip = fabs( Tbvjnin - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == NoK 

data [ i 3 [0] = qb; data [ i ] [ 1 ] = Tb; data [ i ] [2] = 1; 
data [ i ] [3] = 0; data [i] [4] = C; data [i] [5] = 1; 

Ta = Tb; qa = qb; ++i;> 



a = bl * length; 

while (i < = N && overflow == NoK 
segment jiiff = 10; 

a_width = edge_width * a / ( bl - length); 
b_width = edge_width * ( a + db ) / ( bl - length); 

while ( fabs( segment jji ff ) > ( Tb * error) / N && overflow == NoK 

ua = u(a) ; ub = u(a + db); 

A = A_( ) ; B = B_( ) ; C = C_(); D = D_(); 
temp = A*Ta + B* qa; 
qb = C*Ta + D* qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment jii f f ; 

if (Tb > Tbase + 50 | | Tb < OK 

overflow = Yes; possi ble_overf low = Yes; 
if (Ttip < Tbv_min ) Tbvjnin = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == NoK 

data [i] [0] = qb; data [i] [1] = Tb; data [i] [2] = A; 
data [ i ] [33 = B; data C i ] [4] = C; data C i ] [5] = D; 

Ta = Tb; qa = qb; a = a + db; ++i;}} 



if (number_of_tri es > 29K 

clrscrO; cprintf ("Fatal Error: Unable to Compute"); 

getch(); exit(0);} 

if (overflow == NoK 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overf low == Yes){ 
inc = inc / 2; 

if (overal l_di f f >0) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overal l_di ff ; 
number_of_tries++; 
gotoxy(27 / 24); 

cprintf("N = %d Ta = %4.1f Tb = %4.1f ",number_of_tries,w_T(T,Ttip),w_T(T # Tk)) 
if (Ttip < Tave) Ttip = Tave;}} 
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for (i = N - 1; i >= 1; i - - K 

Al = A * data [i ] [2] ♦ B * data[i][4]; B1 = A * data[i] [3] + B * data[i][5); 

Cl = C * data [i ) [2] + D * data[i][4]; D1 = C * data[i] [3] + 0 * data Ci 3 C53 ; 

A = Al; B = B1 ; C=C1; D=D1;> 

length = length + db; 

if (S == 1){ 

out = fopen( f i le_name, "w"); 
for (i = 1; i <= N; i++) 

fprintf (out , "°/ 0 e, %e, °/oe\n'\ w_L(d, i * length / N ) , w T(T,data [i] [1] ) , w H(H, 

data[i] CO] )); 

fclose(out);} 

Ta = Ttip; 

Ap = .5 * length * base_width; 
mass = Ap * depth * density * 1000; 

PN = 2 * sigma * epsilon * pow( length, 2) * pow(Tbase,3) / ( k*base_wi dth ) ; 

eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase, 4 ) ) ; 

result s( ); 

exit(0); 

return;} 

j ★ ★ ★★★★★★****★★★*★★★******★*★★*★★****•**★★**************★•*★***★★★*★**★**** 

OUTPUT PROCEDURE 

★ ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★A***** + j 



void resul ts(void)C 



c l rscr( ) ; 

cprintf ( "\a THE LONGITUDINAL FIN OF TRIANGULAR PROF I LE\r\n" ) ; 

cpr i nt f ( " ==================================================================\r\n" ) ; 

cprintfC Density = %4.2f kg / dm~3\r\n" ,densi ty) ; 

cprintf(" Fin Depth = 11 ); p_L(L, depth); 

cprintfC" Total Mass = "); p_M(M, mass); 

cprintfC Fin Length = "); p_L(L, length); 

cprintfC" Emissivity = %4.2f\r\n", epsilon); 

cprintfC Profile Area = "); p_A(X,Ap); 

cprintfC Profile Number = %4 .3f \r\n", PN) ; 

cprintfC Base Heat Flow = "); p_H(H, qb); 

cprintfC Fin Efficiency = %4 .3f \r\n", eta) ; 

cprintfC Fin Base Width = "); p_L(d, base_width); 

cprintfC Base Temperature = 11 ); p_T(T,Tb); 

cprintfC Edge Temperature = "); p_T(T,Ta); 

cprintfC Thermal Conductivity = "); p_K(K, k); 

cprintfC Environmental Parameter = %4.2f\r\n", V); 

cprintfC ==================================================================\ r \n"); 

cprintfC Transmission Parameter Matrix \r\n\r\n"); 

cprintfC A = %4.4f B = ", A); p_Z(Z,B); 

cprintfC C = "); p_Y(Y,C); 

cprintfC D = %4.4f", D); 

getch( ); 

return;} 



jit *★★★★***★★★★★***★★★★★★**★★★★*******★★*★★★★★***★★*★****★★★★★**★*★★★★★★★★ 

INPUT PROCEDURE 

A******************************************************************'******* it J 



void input (void)C 

clrscr( ); _set cursor type(_NORMALCURSOR ) ; 



cprintf ("\r\n THE LONGITUDINAL FIN OF TRIANGULAR PROF I LE\r\n") ; 

cprintf ("\r\n ==================================================================\r\n") ; 

cprintfC Fin Depth "); depth = i_L(L); 

cprintfC Fin Length "); length = i_L(L); 



87 



cprintf(" Fin Base Width "); base_width = i_L(d); 

cprintfC" Emissivity = "); scant ("%f", Sepsilon); 

cprintfC" Density ( kg / dnT3 ) = "); scant ("%f n ( ^density); 

cprintt(" Base Temperature "); Tbase = i_T(T); 

cprintt(" Thermal Conductivity "); k = i K ( K ) ; 

cprintTC" Environmental Parameter (0 to .8) = "); scant("%t", &V); 

it (S == 1) C 

cprintt(" Output Filename is = "); 

scantC'Xs", tile_name);> 

Cpri ntt C"\r\n =============-===================================================="); 



it ( epsilon > 1.0 || epsilon <= 0){ 

cprintt("\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 

getch(); exit(0);) 

it ( V < 0 j | V > 0.8K 

cprintt("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 

getchO; exit(0);) 

edge_width = base_width / 1e2; 

_setcursortype(_NOCURSOR ); 
return; } 

j * ★★★★****★★★★★*★*****★****★***★********★★★********★★★**★***★★***★★**★*** 

GLOBAL FUNCTIONS 

************************************************************************** ★ j 



double hr(void){ 

tloat Tav = (Tb + Ta) / 2; 

return(sigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)));) 
double m(void)C 

return( hr() / (k * sin(tin_taper)));> 

double u(double x){ 

return( 2 * sqrt(x * m( ) ) );) 

double A_(void){ 

return( ua * ( 1 1 (ua)*K0(ub) + I0(ub)*K1 (ua) ) );} 

double B_(void){ 

return ( ub*ub*( I0(ub)*K0(ua)- 10(ua)*K0(ub)) / ( 2*k*b_width*depth*m( ) ) );) 
double C_(void){ 

return ( 2*k*a_width*depth*m( )*ub*( 1 1 (ub)*Kl (ua) - I 1(ua)*Kl (ub) )/(ua) );) 
double D_(void)( 

return( ub*( I0(ua)*Kl (ub)+ I 1(ub)*K0(ua) ) );) 
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j* A************************************************************************** 

PROFILE AREA Analysis Program 
Longitudinal Fin of Rectangular Profile 

ft***************************************************'*'*'*’*'*’***'**'************** ★ j 



U include <stdio.h> 
tfinclude <stdlib.h> 

# include <math.h> 

# include <conio.h> 

U include "units.h" 

^define Yes 1 

#define No 0 

^define N 100 

tfdefine sigma 5.66961e-8 

char f i le_name [20] ; 

double length, width, depth, start i ng_width , endi ng_width, width_inc; 

double A, B, C, D, Ta, Tb, qb=0; 

float V, K2, epsilon, k, Tbase, density, mass, Ap, error, eta; 

int T, L, M, S, H, d, K, i, X, Z, Y; 

void inputO, computeO; 

j* *********************************************************************** 



MAIN PROGRAM 

************************H*************************************iii**i****** * j 



int mainOC 

FILE *data_in, *out; 

data_in = fopenC'uni ts.ext","r"); 

f scanf (data_in, "%d %d %d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 
f cl ose(data_in) ; 
error = error / 100; 

inputO; 

width_inc = (ending_width - start ing_width ) / 100; 
width = start ing_width; 

while (width < endi ng_width )C 

compute ( ); 

while (qb == OK 

width = 1 .001 * width; 
compute ( );) 

gotoxyC 19, 23) ; 

cprintf ("width = %4.1f length = %4.1f heat = %4.1f", w_L(d, width) , w_L(L , length) , 
w_H ( H , qb)); 

if (S == 1)C 

out = f open( f i le_name, "a"); 

fprintf(out, "%e, %e, %e, %e, %e, %e, %e, %e, %e, %e\n", w_L(d, width) , 

w_T ( T , Tb) , w_T ( T , T a ) , eta, w_H(H, qb), w_M(M,mass), A, w_Z(Z,B), w_Y(Y,C), D) 
fclose(out );} 

width = width + width_inc;) 



gotoxy(21 ,24); cl reol( ); cprintf ( "\aCompleted \r\n"); 

clreolO; 

getch( ); 
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exit(O); 
return 0;} 

jk kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

COMPUTE ROUTINE 

************************************************************************** k j 



void compute(void){ 

double A_( ) , B_(), C _(), D_( ) , m(), hr(), Yb(); 

double inc f overalljdiff = 100, Tmax = 0, Tbv_min = 100000, data [N+1] [6] ; 

double Al , B1, Cl, 01, Ttip, temp, qa, Tave, segment_di f f ; 

int overflow, possible_overf low = No, number_of_tries=0; 

length = Ap / width; 

K2 = 2 * epsilon * sigma * V * pow(Tbase, 4) ; 

Tave = Tbase * pow (V, 0.25) + 0.0001; 

if (V==0 ) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2) 

else 

Tt ip = .0001 + Tave; 
inc = .5 * Ttip; i = 1; 

while ( fabsCoveral l_di f f ) > ( Tbase * error) || i <= N K 

overflow = No; i =1; Ta = Ttip; Tb = Ta; 

qa = sigma * width * epsilon * pow(Ta, 4) * depth; 

while (i <= N && overflow == NoH 

segment_diff = 10; 

while ( fabs(segment_di f f ) > ( Tb * error) / N && overflow == No){ 

A = A_(); B = B_( ) ; C = C_(); D = D_(); 
temp =A*Ta+B*qa; 
qb =C*Ta+D* qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_di f f ; 

if (Tb > Tbase + 50 J J Tb < OK 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;)) 

if (overflow == No)( 

data [i 3 [03 = qb; data [i] [13 = Tb; data[i)[2] = A; 
dataCi] [3] = B; data [i 3 [43 = C; data[i] [5] = D; 

Ta = Tb; qa = qb; ++i;)) 

if (number_of_t ri es > 29){ 
qb=0; return;) 

if (overflow == NoH 

if (Tmax < Ttip) 

if (Tb < Tbase ) Tmax = Ttip; 

overalljdiff = (Tb - Tbase) / 2; 

if ( possible_overf low == YesK 
inc = inc / 2; 

if (overalljdiff >0) Ttip = Ttip - inc; 
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else 



else Ttip = Ttip + inc;> 



Ttip = Ttip - overal l_di f f ; 
number_of_tri es++; 
gotoxy(27, 24) ; 

cprintf("N = %d Ta = %4.1f Tb = %4 . 1 f ",number_of_t ri es, w_T ( T , Tt ip) , w_T ( T , Tb) ) 
if (Ttip < Tave) Ttip = Tave;}} 

for ( i = N - 1; i >= 1; i - - H 

A1 = A * datati] [2] + B * datati] [4] ; B1 = A * datati] [3] + B * datali] [5] ; 

Cl = C * dataCi] [2] + D * datati] [4]; D1 = C * datati] [3] + D * datati] [5] ; 

A = A 1 ; B = B 1 ; C=C1; D=D1;} 

Ta = Ttip; 

eta = qb / ( 2 * sigma * epsilon * length * depth * powCTbase,4) ); 

mass = width * length * depth * density * 1000; 

return;} 

Jk ★★★★A****************************************************************** 

INPUT PROCEDURE 

************************************************************************** * I 



void input (voidK 

cl rscr( ) ; _setcursortype(_NORMALCUR$OR ) ; 



cprintf ("\r\n THE LONGITUDINAL FIN OF RECTANGULAR PROFILE"); 

cpr intf ( "\r\n ==================================================================\r\n" ) ; 

cprintfC" Fin Depth "); depth = i_L(L); 

cprintf(" Starting Fin Width "); start ing_width = i_L(d); 

cprintfC" Ending Fin Width "); ending_width = i_L(d); 

cprintfC" Profile Area "); Ap = i_A(X); 

cprintfC" Emissivity = scanf("%f", Sepsilon); 

cprintfC" Density C kg / dm A 3 ) = "); scanf("%f", ^density); 

cprintfC" Base Temperature "); Tbase = i __T C T ) ; 

cprintfC" Thermal Conductivity "); k = iJCCK); 

cprintfC" Environmental Parameter CO to .8) = "); scanfC"%f", &V); 

if (S == 1) C 

cprintfC" Output Filename is = "); scanf("%s", file_name);} 

cprintf C "\n ==================================================================\ r \n") ; 



if C epsilon > 1.0 j| epsilon <= OK 

cpr intf ( "\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 

getchC); exit(0);} 

if C V < 0 | i V > 0.8)0 

cpr intf C"\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 

getchC); exit CO);} 

_setcursortype(_NOCURSOR) ; 
return;} 

jk kkk*k*kkkkkkkkk*kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk*kkkkkkkkkkkkk 

GLOBAL FUNCTIONS 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk k j 



double hr(void){ 

float Tav = C T b + Ta) / 2; 

returnCsigma * epsilon * pow(Tav, 3) - CK2 / (2 * Tav)));} 
double m( voidK 

returnC sqrtC C 2.0 * hrO) / C k * width) ) );} 
double Yb(voidK 
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return( sqrt( 2.0 * hr() * k * width) * depth );) 



double A_(voidK 

returnC cosh( m( ) * length / N 

double B_(voidX 

return ( sinh( m( ) * length / N 

double C_(voidX 

return ( sinh( m( ) * length / N 

double D_(void){ 

returnC cosh( m( ) * length / N 



);> 

) / Yb() );> 
) * Yb() );} 
);> 
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y* ************************************************* 



************************ 



PROFILE AREA Analysis Program 
Longitudinal Fin of Trapezoidal Profile 

**************************************************************************** * j 



#include <stdio.h> 
#include <stdl ib.h> 
# include <math.h> 

# include <conio.h> 
^include "units. h" 
tfinclude "bessel ,h" 



#def i ne 
#def ine 
#def ine 
#def ine 



Yes 1 

No 0 

N 100 

sigma 5. 66961 e- 8 



char 

double 

double 

float 

int 

void 



f i le_name [20] ; 

length, edge_width, base_width, a_width, b_width, depth, start i ng_width, er>ding_width 
A , B, C, D, Ta, Tb, qb=0, ua, ub, fin_taper, a, width__inc; 

V, K2, epsilon, k, Tbase, taper_ratio, density, mass, Ap, error, eta; 

T, L, M, S, H, d, K, i, X, 2, Y; 
inputC), computeC); 



y* *********************************************************************** 



MAIN PROGRAM 



************************************************************************** *y 



int main(){ 

FILE *data_in, *out; 

data_in = fopenC'uni ts.ext", "r"); 

fscanf (data_in,"%d %d %d %d %d %d %d %d %d %d %f", &T,&L,&M,&S,&H, &d,&K,&X,&Z,&Y, terror); 

fclose(data_in); 

error = error / 100; 

i nput ( ) ; 

width_inc = (ending_width - start ing_width ) / 100; 
base_width = start ing_width; 

while (base_width < ending_width){ 

computeC ); 

while (qb == 0)( 

base_width = 1.001 * base_width; 
computeC );} 

gotoxyC 19, 23) ; 

cprintf ("width = %4.1f length = %4.1f heat = %4.1f", w_L(d, base_width) , w_L(L, length) , 
w_H( H , qb)); 

if (S == 1)( 

out = fopen(f i le_name, "a"); 

fprintfCout, "%e, %e, %e, %e, %e, %e, %e, %e, %e, %e\n", w L(d,base_width), 
w_T ( T , Tb) , w_T ( T , T a ) , eta, w_H(H, qb), w_M(M,mass), A, w_Z(Z,“B), wJT(Y,C), D); 
fclose(out);} 

base_width = base_width + width_inc;> 



gotoxyC 21 , 24);clreol ( ); cprintf ("\aCompleted \r\n"); 
cl reol ( ) ; 
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getchC ) ; 
ex i t ( 0 ) ; 
return 0;> 

I* A********************************************************************** 

COMPUTE ROUTINE 

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★***************************** * i 



void compute( voidK 



double 
double 
double 
double 
i nt 



A_(), B_( ) , C_(), D_( ) , m( ) , hr(), u(double x); 
lengthy inc, bl , db; 

inc, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 100000, data [N+1] [6] ; 
A1 , Bl, Cl, D1 , Ttip, temp, qa, segment_di f f ; 
overflow, poss i ble_overf low = No, number_of_tri es=0; 



edge_width = base_width * taper_ratio; 
length = 2 * Ap / ( base_width + edge_width ); 
bl = length / ( 1 - taper_ratio ); 
fin_taper = atan2( base_width, 2*b1 ); 
db = length / N; 



K2 = 2 * epsilon * sigma * V * pow(Tbase,4 ) ; 
Tave = Tbase * powCV, 0.25) + 0.0001; 



if (V==0) 

Ttip = .8625 * Tbase * sqrtC 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = Tave + 0.0001 ; 
inc = .5 * Ttip; i=1; 

while (fabs(overal l_diff ) > ( Tbase * error) ]] i <= N K 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * edge_width * epsilon * pow(Ta, 4) * depth; 

a = bl - length; 



while (i <= N && overflow == NoK 



segment_diff = 10; 

a_width = edge_width * a / ( bl - length); 
b_width = edge_width *(a+db)/(b1 - length); 

while ( f abs(segment_di f f ) > ( Tb * error) / N && overflow == NoK 



ua = u(a) ; ub = u(a + db); 

A = A_ ( ) ; B = B_( ) ; C = C_(); D = D_(); 
temp =A*Ta+B*qa; 
qb =C*Ta+D* qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segmented if f ; 

if (Tb > Tbase ♦ 50 ] ] Tb < OK 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbv_min ) Tbvjnin = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;)) 

if (overflow == NoK 

data [i ] [03 = qb; datalij [1] = Tb; data [iD [2] = A; 
data [i] [3] = B; data [iD [43 = C; datatij [5] = D; 
Ta = Tb; qa = qb; a = a + db; ++i;)) 



if (number_of_tries > 29 )( 
qb=0; return;) 
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if (overflow == NoH 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overal l_diff = (Tb - Tbase) / 2; 

if ( possible_overf low == Yes)C 
inc = inc / 2; 

if (overal l_di ff > 0 ) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;) 

else 

Ttip = Ttip - overal l_di ff; 
number_of_tri es++; 
gotoxy(27,24); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f ", number_of_t ri es, w_T( T , T t ip) , w_T(T f Tb) ) 

if (Ttip < Tave) Ttip = Tave;}} 

for (i = N - 1; i >= 1; i--)( 

A1 = A * data [i] [2] + B * data [i] [4] ; B1 = A * data[i][3] + B * data[i][5]; 

Cl = C * data[i] [2] ♦ D * data [ i] [4] ; D1 = C * dataCi] [3] + D * data [ i 3 [5] ; 

A = A 1 ; B = B1 ; C=C1; D=D1;> 

Ta = Ttip; 

mass = .5 * ( edge_width + base_width ) * length * depth * density * 1000; 
eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase, 4)); 
return;) 

I ★ ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★A'*** 

INPUT PROCEDURE 

★ ★a************************************** ********************************* * j 



void input(void)( 

clrscr( ); _setcursortype(_NORMALCURSOR ) ; 



cprintf (" THE LONGITUDINAL FIN OF TRAPEZOIDAL PROFILE"); 

cprint f ( "\r\n ==================================================================\r\n"); 

cprintf(" Fin Depth "); depth = i_L(L); 

cprintf(" Starting Base Fin Width "); start ing_width = i_L(d); 

cprintf (" Ending Base Fin Width "); ending_width = i_L(d); 

cprintf(" Taper Ratio = "); scanf("%f", &taper_rat io); 

cprintfC Profile Area "); Ap = i_A(X); 

cprintf(" Emissivity = "); scanf ("%f ", Sepsilon); 

cprintfC Density ( kg / dm A 3 ) = "); scanf C^f", &density); 

cprintfC Base Temperature "); Tbase = i_T(T); 

cprintfC Thermal Conductivity "); k = i_K(K); 

cprintfC Environmental Parameter (0 to .8) = "); scanf ("%f" # &V); 

if (S == 1) C 

cprintfC Output Filename is = "); scanf ("%s", file_name);) 

cprintf ("\n ==================================================================\r\ n »); 



if ( taper_ratio > .99 |] taper_ratio <= 0.0) ( 

cprintf ("\r\n\r\nFATAL ERROR: Taper Ratio is not between 0 and 1"); 

getchO; exit(0);) 

if ( start ing_width >= ending_wi dth)( 

cprintf ("\r\n\r\nFATAL ERROR: Starting Base Fin Width >= Ending Base Fin Width"); 

getchO; exit(0);) 

if ( epsilon > 1.0 ]] epsilon <= OK 

cprintf C\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 

getchO; exit(0);) 

if ( V < 0 | | V > 0.8K 

cprintf ("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 

getchO; exit(0);) 
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_setcursortypeC_NOCUR$OR ); 
return; ) 

I* ★★*★*★**★★★★★★★★★★*★****★★***★*★*★*****★***********★*****★**★★★**★***** 

GLOBAL FUNCTIONS 

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★A * j 



double hrCvoidX 

float Tav = (Tb + Ta) / 2; 

returnCsigma * epsilon * powCTav, 3) - (K2 / (2 * Tav)));) 
double mCvoidX 

returnC hr() / (k * sin(f in_taper) ));) 

double uCdouble x){ 

returnC 2 * sqrt(x * m( ) ) );) 

double A_(void)( 

returnC ua * C 1 1 Cua)*K0Cub) + I0Cub)*Kl (ua) ) );} 

double B_(void){ 

return C ub*ub*C 1 0Cub)*K0Cua) - 1 0Cua)*K0Cub) ) / C 2*k*b_width*depth*m( ) ) );) 
double C_(voidX 

return C 2*k*a_width*depth*mC )*ub*C 1 1 Cub)*Kl Cua) * 1 1 Cua)*K1 Cub) )/Cua) );) 
double D_Cvoid){ 

returnC ub*C 1 0Cua)*K1 Cub)+ 1 1 Cub)*K0Cua) ) );) 



96 



j * *************************************************************************** 



PROFILE AREA Analysis Program 
Longitudinal Fin of Triangular Profile 

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★a ★ j 



^include <stdio.h> 
^include <stdlib.h> 
# include <math.h> 
tfinclude <conio.h> 
tfinclude "units. h" 
#include "bessel.h" 



#def ine 
#def ine 
#def ine 
#def ine 



Yes 1 

No 0 

N 99 

sigma 5. 66961 e- 8 



char f i le_name [20] ; 

double length, width_edge, base_width, a_width, b_width, depth, start ing_width , ending_width 

double A, B, C, D, Ta, Tb, qb=0, ua, ub, fin_taper, a, inc_width; 

float V, K2, epsilon, k, Tbase, density, mass, Ap, error, eta; 

int T, L, M, S, H, d, K, i, X, Z, Y; 

void input(), computeO; 



j-k ****★***********★★****★★★★******'*★*'***★*★**★*★**★★**★★★★*★★*★★★★★★★★★★★ 



MAIN PROGRAM 



★★*★*★**★****★*** ********************** A********************************** ★ j 



int main(K 

FILE *data_in, *out; 

data_in = fopen("units.ext","r"); 

f scanf (data_in, "%d %d % d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error = error / 100; 

input ( ); 

inc_width = (ending_width - start ing_width ) / 100; 
base_width = start ing_width; 

while (base_width < ending_width){ 

compute ( ); 

while (qb == 0){ 

base_width = 1.001 * base_width; 
compute( ) ; > 

gotoxy( 19,23); cprintfC "width = %4.1f length = %A.1f heat = %4.1f", w_L(d,base_width) , 
w_L(L, length), w_H(H, qb)); 

if (S == DC 

out = f open( f i le_name, "a"); 

fprintf(out, "%e, %e, %e, %e, %e, %e, %e, %e, %e, %e\n", w_L(d,base_width) , 

w_T ( T , Tb) , w_T ( T , T a ) , eta, w_H(H, qb), w_M(M,mass), A, w_Z(Z,B), wJT(Y,C), D); 
fc l os e( out ) ;) 

base_width = base_width + inc_width;> 



gotoxy(21 ,24); cl reol ( ) ;cprintf ("\aCompleted \r\n"); 
cl reol ( ) ; 
getch( ) ; 
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exi t(0); 
return 0;> 

I ★ **★★*★*★★★*★★*★★*****★**★****★★★★*★'*★*********★***■*★*★*****★****★★*★★★* 

COMPUTE ROUTINE 



void compute(void)( 

double A_( ) , B_(), C_(), D _(), m(), hr(), u(double x); 

double length_inc, bl , db; 

double inc, Tave # overall_diff = 100, Tmax = 0, Tbvjnin = 100000, data[N+1][6] 

double Al , Bl, Cl, D1, Ttip, temp, qb, segment_di f f ; 

int overflow, possi blejDverf low = No, number_of_tries=0; 

length = 2 * Ap / base_width; 

width_edge = base_width / 100; 

db = length / ( N + 1 ); 

bl = length; 

length = length ■ db; 

fin_taper = atan2( base_width, 2*b1 ); 

K2 = 2 * epsilon * sigma * V * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2) 

else 

Ttip = Tave + 0.0001 ; 
inc = .5 * Ttip; i=1; 

while (fabs(overal l_dif f ) > ( Tbase * error) ] ] i <= N )( 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 
segment_diff = 10; 

while (f abs(segment_di f f ) > ( Tb * error) / N && overflow == No)( 

ub = u(db); 

temp = Ta / I0(ub); 

C = 2 * k * width_edge * depth * pow(n(),2) * 1 1 ( ub ) / (ub * I0(ub)) 
qb = C * Tb; 

segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_di f f ; 

if (Tb > Tbase +50 [ | Tb < 0)( 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbv_mi n ) Tbv_mi n = T t i p; 
inc = (Tbvjnin - Tmax) / 2; 

Ttip = fabs( Tbvjnin - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == No){ 

data [i] [0] = qb; da t a [ i 3 [ 1 } = Tb; dataCi) [2] = 1; 
data C i 3 [33 = 0; data Ci3 [43 = C; data [ 13 [53 = 1; 
a = bl - length; 

Ta = Tb; qa = qb; ++i;} 

while (i <= N && overflow == NoK 

segment_diff = 10; 

a_width = width_edge * a / ( bl - length); 
b_width = width_edge * ( a + db ) / ( bl - length); 

while (fabs(segment_di f f ) > ( Tb * error) / N && overflow == NoK 
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ua = u(a) ; ub = u(a + db); 

A = A_( ); B = B_( ); C = C_(); D = D_(); 
temp =A*Ta+B* qa; 
qb = C * Ta + D * qa; 
segmented! ff = ( temp • Tb) / 2; 

Tb = temp + segment_di f f ; 

if (Tb > Tbase -*-50 j | Tb < OX 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == NoX 

data [i ) [0] = qb; data [ i ] [1 ] = Tb; data [i D [2] = A; 
data [i ) [3] = B; data [ i ) [4] = C; data [i ] [5] = D; 

Ta = Tb; qa = qb; a = a + db; ++i;>> 

if (number_of_tr i es > 29X 
qb=0; return;} 

if (overflow == NoX 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possi ble_overf low == YesX 
inc = inc / 2; 

if (overatl_diff > 0 ) Ttip = Ttip - inc; 
else Tt ip = Ttip + inc;} 

else 

Ttip = Ttip - overal l_di f f ; 
number_of _t r i es++ ; 
gotoxy(27,2A); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f ", number_of_tr i es, w_T ( T , T t ip) , w_T ( T , Tb) ) 

if (Ttip < Tave) Ttip = Tave;}} 

for (i = N - 1; i >= 1; i-- X 

A1 = A * dataCi] [2] + B * dataCi] [4]; B1 = A * data[i] [3] + B * data[i] [5] ; 

Cl = C * dataCi] [2] + D * dataCi] [4]; D1 = C * data [ i ] C3] + D * data[i][5]; 

A = Al; B = Bl; C=C1; D=D1;} 

length = length + db; 

Ta = Ttip; 

mass = .5 * base_width * length * depth * density * 1000; 

eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase,4) ); 

return;} 

!+ *★★**★★**★**★★**★**★★★★***★★★*★★★★★★★★★★*★★★★★★★★********★******★*★**** 

INPUT PROCEDURE 

***************************************** ******************* ************** * j 



void input(void){ 

clrscr( ); _set cursor type(_NORMALCURSOR); 



cprintfC THE LONGITUDINAL FIN OF TRIANGULAR PROF I LE\r\n" ) ; 

cprintf (" =========================rr=======================================\ r \n"); 

cprintfC Fin Depth "); depth = i_L(L); 

cprintfC Starting Base Fin Width "); start i ng^width = i_L(d); 

cprintfC Ending Base Fin Width "); ending_width = i_L(d); 

cprintfC* Profile Area "); Ap = i _ A ( X ) ; 

cprintfC Emissivity = "); scanf("%f", Sepsilon); 

cprintfC Density ( kg / dm A 3 ) = *'); scanf("%f", Sdensity); 

cprintfC Base Temperature "); Tbase = i_T(T); 

cprintfC* Thermal Conductivity *'); k = i_K(K); 
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cprintf ( 



Environmental Parameter (0 to .8) = "); scanf("%f" # &V); 



cprintfC Output Filename is = "); scanf("%s", file_name);> 

cprintf ( "\n =rs=s===r===========r=================rr==rr=rsssssrss=====s======\ r \ n "); 



if ( start ing_width >= ending_width){ 

cprintf ("\r\n\r\nFATAL ERROR: Starting Base Fin Width >= Ending Base Fin Width") 

getch(); exit(0);} 

if ( epsilon > 1.0 jj epsilon <= OK 

cprintf ("\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 

getch(); exit(0);} 

if ( V < 0 | | V > 0.8K 

cprintf ("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8") 

getch(); exit(0);) 

_setcursortype(_NOCURSOR ) ; 
return;} 

I* **★********★*★*********★**★******★*★★★★******★★★★**********★*★**★**★*** 

GLOBAL FUNCTIONS 

★A************************************************************************ ★ j 



double hr(void){ 

float Tav = (Tb + Ta) / 2; 

return(sigma * epsilon * pow(Tav, 3) - (K2 / ( 2 * Tav)));} 
double m(void){ 

return( hr() / (k * sin(f in_taper) ) );} 

double u(double x)( 

return( 2 * sqrt(x * m( ) ) );} 

double A_(void){ 

return( ua * ( 1 1 (ua)*K0(ub) + 1 0(ub)*K1 (ua) ) );} 

double B_( voidX 

return ( ub*ub*( I0(ub)*K0(ua)-l0(ua)*K0(ub)) / ( 2*k*b_width*depth*m( ) ));} 
double C_(voidK 

return ( 2*k*a_width*depth*m( )*ub*( 1 1 (ub)*K1 (ua) - 1 1 (ua)*K1 (ub) )/(ua) );} 
double D_(void){ 

return( ub*( I 0(ua)*K1 (ub)+ 1 1 (ub)*K0(ua) ) );} 

jit ★***★★***★*★*★★****★***★*******★*★************★*★******★*★★*★★★★★**■*★*★★*★* 



OPTIMIZE Analysis Program 
Longitudinal Fin of Rectangular Profile 



ft****************************************************************'*'********** ★ j 



# include 


<stdio.h> 




# include 


<stdl ib.h> 




# include 


<math .h> 




#include 


<conio.h> 




# include 


"uni ts.h" 




#def i ne 


Yes 


1 


#def ine 


NO 


0 


#def i ne 


N 


100 


#def ine 


s i gma 


5.66961e-8 


double 


length 


, width, depth, 



100 



float 
i nt 
void 



Tbase, heat, PN, eta, V, K2, epsilon, k, density, mass, Ap, error; 
T, L, M, S, H, d, K, i, X, 2, Y; 
inputO, computeO, resultsC); 



jit *********************************************************************** 



MAIN PROGRAM 

★ ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★A** * j 



int main()( 

FILE *data_in; 

data_in = f openC'uni ts.ext", M r M ); 

f scanf (data_in, "%d %d %d %d %d %d %d %d %d %d %f",S,T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error = error / 100; 

inputO; 

heat = heat / ( 1 - V); 

Ap = (1.6178 * pow(heat/depth , 3)) / ( k * pow(epsi Ion, 2) * pow(Tbase,9) * pow(sigma,2) ) 
width = (1.8648* pow(heat/depth , 2)) / ( k * sigma*epsi lon*pow(Tbase, 5) ); 

computeO; 
resul ts( ); 
ex i t ( 0 ) ; 
return 0;> 

j-k 

COMPUTE ROUTINE 

★************************************************************************* * j 



void compute(voidK 

double A_( ) , B_( ) , C_( ) , D_( ) , m(), hr(), Yb(); 

double inc, overall_diff = 100, Tmax = 0, Tbv_min = 100000, data [N+1] [6] ; 

double A1 , B1, Cl, 01, Ttip, Tave, temp, qa, segmentjdi f f ; 
int overflow, possi ble_overf low = No, number_of_tri es=0; 

length = Ap / width; 

K2 = 2 * epsilon * sigma * V * pow(Tbase,4 ) ; 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if ( V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = .0001 + Tave; 
inc = .5 * Ttip; i = 1; 

while (fabs(overal l_di f f ) > ( Tbase * error ) J | i <= N ) C 

overflow = No; i = 1; Ta = Ttip; Tb - Ta; 

qa = sigma * width * epsilon * pow(Ta, 4) * depth; 

while (i <= N && overflow == NoX 

segment _diff = 10; 

while (f abs(segment_di f f ) > ( Tb * error) / N && overflow == No){ 

A = A_( ) ; B = B_( ) ; C = C_(); D = D_(); 
temp =A*Ta+B*qa; 
qb = C * Ta + D * qa; 
segment_diff = ( temp - Tb) / 2; 
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Tb = temp + segmented i f f ; 

if (Tb > Tbase + 50 } J Tb < OK 

overflow = Yes; possible_overf low = Yes; 
if ( T t i p < T bv_mi n ) T bv_mi n = Ttip; 
inc = (Tbv_min - Tmax) / 2 ; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;)} 

if (overflow == NoX 

data [i] [0] = qb; data [i ] [ID = Tb; data [ i D [2] = A; 
data C i ] [3] = B; datali] [4] = C; data [ i ] [5] = 0; 

Ta = Tb; qa = qb; ++ i ;}} 

if (number_of_tri es > 29 X 

clrscrO; cpr i ntf ("Fatal Error: Unable to Compute"); 

getch(); exit(0);} 

if (overflow == No)( 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overf low == YesX 
inc = inc / 2; 

if (overal l_di f f > 0 ) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overal l _di ff ; 
number_of _t r i es++; 
gotoxy(27,24); 

cprintfCN = %d Ta = %4.1f Tb = %4.1f " f number_of_t ries, w_T(T , Tt i p) , w_T (T , Tb) ) 

if (Ttip < Tave) Ttip = Tave;}} 

for (i = N - 1; i >= 1; i--X 

A1 = A * data [ i ] [2] + B * datali] [4]; B1 = A * dataCi] [3] + B * data[i] [5] ; 

Cl = C * datali] [2] + D * datali] [4]; D1 = C * datali] [3] + D * datali] [5] ; 

A = Al; B = B1 ; C=C1; 0=01 ; } 

Ta = Ttip; 

eta = qb / (2 * sigma * epsilon * length * depth * pow(Tbase,4)); 

PN = length * sqrt(sigma * epsilon * pow(Tbase,3) / (k * width)); 
mass = width * length * depth * density * 1000; 
return; } 

j-k ★★****★★★★★★**★*★★★***★********★*'*★★★★★★★*★****★★★★★★★★★*★★★★★★★★★★★★★★ 

OUTPUT PROCEDURE 

a ********************************************************************** it j 



void results(void){ 



c l rscr( ) ; 

cprintf ("\a OPTIMUM PARAMETERS"); 

cprintf ("\r\n THE LONGITUDINAL FIN OF RECTANGULAR PROFILE"); 

cprintf ("\r\n ====================================:==============================\r\n") ; 

cprintf (" Density = %4.2f kg / dm A 3\r\n", densi ty); 

cprintf (" Fin Depth = "); p_L(L, depth); 

cprintfC Total Mass = "); p_M(M # mass); 

cprintf (" Fin Height = "); p_L(L, length); 

cprintfC 1 Fin Width = "); p_L(d, width); 

cprintfC Emissivity = %4.2f\r\n", epsilon); 

cprintfC Profile Area = "); p_A(X, Ap); 

cprintfC Profile Number = %4 .3f \r\n", PN ) ; 

cprintfC Base Heat Flow = "); p_H(H, qb); 

cprintfC Fin Efficiency = %4 .3f \r\n", eta); 

cprintfC Base Temperature = "); p_T(T,Tb); 
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cprintfC Edge Temperature = “); p_T(T,Ta); 

cprintfC" Thermal Conductivity = "); p_K(K, k); 

cprintfC Environmental Parameter = %4.2f\r\n", V); 

cprintfC 11 ==================================================================:\ r \n" ) ; 

cprintfC Transmission Parameter Matrix \r\n\r\n"); 

cprintfC A = %4.4f B = ", A); p_Z(2,B); 

cprintfC" C = "); p_YCY,C); 

cprintfC* D = %4.4f", D); 

getchC ) ; 

return; } 



jk kk*kkkkk*kkkkk*kkkkkkkkkkkkkkkkkkkkkk*kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

INPUT PROCEDURE 



*★*★★*★★★★★*★★★★★***★*★★**★*★★*★****★★*★★*★**★**★**★*★★****★*★*★★★**★★★★★★ ★ j 



void inputCvoid)( 

clrscrC ); _setcursortypeC_NORMALCURSOR ); 



cprintf C\r\n THE LONGITUDINAL FIN OF RECTANGULAR PROF ILE\r\n") ; 

cpri ntf C\r\n ------------------------------------------------------------ ------\ v \^) 

cprintfC Fin Depth "); depth = i_LCL); 

cprintfC" Heat Dissipation "); heat = i_HCH); 

cprintfC Emissivity = "); scanfC%f", Sepsilon); 

cprintfC" Density C kg / drrT3 ) = "); scanfC"%f", ^density); 

cprintfC Base Temperature "); Tbase = i_T(T); 

cprintfC Thermal Conductivity "); k = i_K(K); 

cprintfC Environmental Parameter CO to .8) = "); scanfC"%f", &V); 

cprint f C\r\n 

if C epsilon > 1.0 j| epsilon <= OK 



cprintf C\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 

getchC ); exit CO);} 

if C V < 0 | | V > 0.8X 

cprintf ("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 

getchC); exitCO);} 

_set cursor t ypeC^NOCURSOR ); 
return; } 

j-k kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk-kkk'kickickicki<ick1< 

GLOBAL FUNCTIONS 

kkkkkkkk*k*k*kk*kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk*kkk***kkkk****** k j 



double hrCvoidX 

float Tav = C T b + Ta) / 2; 

returnCsigma * epsilon * powCTav, 


3) - (K2 / 


(2 * Tav)));} 


double mCvGidX 

returnC sqrtC C 2.0 * hrO) / C k 


* width) ) 


>;} 


double YbCvoidX 

returnC sqrtC 2.0 * hr() * k * width) * depth 


);> 


double A_CvoidX 

returnC coshC mC) * length / N ) 


);> 




double B_(void)( 

return C sinhC mC ) * length / N ) 


/ Yb() );> 




double C_Cvoid)( 

return C sinhC mC) * length / N ) 


* Yb() );> 




double D_Cvoid)( 







103 



return( cosh( m( ) * length / N ) 



);> 
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******** 



y* *************************** 

OPTIMIZE Analysis Program 
Longitudinal Fin of Trapezoidal Profile 
************************************************** 



****************** 



***************** 



# include 


<stdi o.h> 




# include 


<stdl ib.h> 




#i nc lude 


<math ,h> 




# include 


<coni o.h> 




#i nclude 


"uni ts ,h" 




# i nc l ude 


"bessel .h" 




#def i ne 


Yes 


1 


#def i ne 


No 


0 


#def i ne 


N 


100 


#def i ne 


s i gma 


5 .66961 e-8 


char 


file name [20] ; 


double 


l ength 


, edge_width 


double 


A, B, 


C , D , T a , a , 


float 


K2, taper_ratio, V 


int 


T, L, 


M, S, H, d, 1 


void 


inputO, computeO 


FILE 


*data_ 


in, *out; 



*********** 



*********** 



******************** 



MAIN PROGRAM 



******************************************************************* 



int main(H 

data_in = f open( "uni ts . ext" , "r" ) ; 

fscanf (data_in,"%d %d %d %d %d %d %d %d %d %d %f",&T ,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 
f close(data_i n) ; 
error = error / 100; 

input C ) ; 

heat = heat / (1 - V); 

Constant_3 = 0.3881*pow(taper_rat io,2) - 1 . 0324*taper_rat i o + 2.7522; 

Ap = ((0.5147 * taper_ratio + 1.1382) * pow(heat/depth , 3)) / ( k * pow(epsi Ion, 2) * 

pow(Tbase,9) * pow(s i gma , 2) ) ; 

base_width = (Constant_3 * pow(heat/depth , 2)) / ( k * sigma*epsi lon*pow(Tbase,5)); 
edge_width = taper_ratio * base_width; 

compute ( ); 
results( ); 
exit(0); 
return 0;> 

!* *********************************************************************** 

COMPUTE ROUTINE 

************************************************************************** * j 

void compute ( voidX 

double A_( ) , B_(), C_(), D_(), m() ( hr(), u(double x); 

double length_inc, bT, db; 

double inc, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 100000, data [N+1] [6] ; 

double A1 , B1, Cl, D1, Ttip, temp, qa, segmented! ff; 

int overflow, possible_overf low = No, number_of_tr i es=0; 
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Length = 2 * Ap / ( base_width + edge_width ); 
bl = length / ( 1 - taper_ratio ); 
fin_taper = atan2( base_width, 2*b1 ); 
db = length / N; 

K2 = 2 * epsilon * sigma * V * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if ( V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = Tave + 0.0001; 
inc = .5 * Ttip; i=1; 

while ( f abs(overal l_di f f ) > ( Tbase * error) J ] i <= N H 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * edge_width * epsilon * pow(Ta, 4) * depth; 

a = bl - length; 

while (i <= N && overflow == No)C 
segment_diff = 10; 

a_width = edge_width * a / ( bl - length); 
b_width = edge_width *(a+db)/(b1- Length); 

while ( fabs(segment_di f f ) > ( Tb * error) / N && overflow == No){ 

ua = u(a) ; ub = u(a + db); 

A = A_( ) ; B = B_(); C = C_(); D = D_(); 
temp =A*Ta+B*qa; 
qb = C * Ta + D * qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_di f f ; 

if (Tb > Tbase + 50 |.| Tb < 0)C 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == NoK 

data C i } CO} = qb; data [ i 3 [ID = Tb; data [ i D [2] = A; 
data [i ] [3] = B; data Cl] C4] = C; data Ci ] [5] = D; 

Ta = Tb; qa = qb; a = a + db; ++i;}} 

if (number_of_tries > 29){ 

clrscr(); cprintf ("Fatal Error: Unable to Compute"); 

getch(); exit(0);} 

if (overflow == No)C 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overf low == YesK 
inc = inc / 2; 

if (overal l_di f f > 0 ) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overa l l_di f f ; 
number_of_tr ies++; 
gotoxy(27, 24 ) ; 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f " , number_of_t r i es, w_T (T , T t i p) , w_T(T , Tb) ) 

if (Ttip < Tave) Ttip = Tave;}} 
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for (i = N - 1; i >= 1; i-*)( 

A1 = A * data C i 3 C2] + B * datati] [4]; B1 = A * datati] [3] + B * datati] t5]; 

Cl = C * data [ i ] [2] + D * datati] [A]; D1 = C * datali] C3] + D * datati] [5]; 

A = Al; B = Bl; C=C1 ; D=D1;> 

if (S == DC 

out = f open( f i L e_name, "w"); 
for (i = 1; i <= N; i++) 

fprintf(out, "%e, %e, %e\n", w_L(d, i * length / N), w_T ( T f data [i] [1] ) , w_H(H, 

datati] [0] ) ); 

f c lose(out );> 

Ta = Ttip; 

Ap = .5 * length * ( edge_width + base_width ); 
mass = Ap * depth * density * 1000; 

PN = 2 * sigma * epsilon * pow( length, 2) * pow(Tbase,3) / (k* base_width); 
eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase,4)); 
return; ) 

I * *********************************************************************** 

OUTPUT PROCEDURE 

**************★★★*******************★****★*******★*****★★**★***★*★★*★*★*■** ★ j 



void resul ts( voi d)( 



c l rscr( ); 

cprintf ("\a\r\n OPTIMUM PARAMETERS"); 

cprintf ("\r\n THE LONGITUDINAL FIN OF TRAPEZOIDAL PROFILE"); 

cpr i ntf ("\r\n ==================================================================\r\n" ); 

cprintf (" Density = %4.2f kg / dm"3\r\n" , densi ty) ; 

cprintfC" Fin Depth = "); p_L(L, depth); 

cprintf(" Total Mass = "); p_M(M, mass); 

cprintfC Fin Length = "); p_L(L, length); 

cprintfC Emissivity = %4.2f\r\n" / epsilon); 

cprintfC Profile Area = "); p_A(X,Ap); 

cprintfC Profile Number = %4.3f\r\n",PN); 

cprintfC Base Heat Flow = "); p_H(H, qb); 

cprintfC Fin Efficiency = %4 .3f \r\n" t eta) ; 

cprintfC Fin Edge Width = "); p_L(d, edge_width); 

cprintfC Fin Base Width = "); p_L(d, base_width); 

cprintfC" Base Temperature = "); p_T(T,Tb); 

cprintfC Edge Temperature = "); p_T(T,Ta); 

cprintfC Thermal Conductivity = "); p_K(K, k); 

cprintfC Environmental Parameter = %4.2f\r\n", V); 

cprintfC =========================r=====r==================================\ r \n") ; 

cprintfC Transmission Parameter Matrix \r\n\r\n"); 

cprintfC A = %A . Af B = ", A); p^Z(Z,B); 

cprintfC C = "); p_Y(Y,C); 

cprintfC D = %4.4f", D); 

getchO; 

return;) 



y* ***★★***********★*********************★*****★****★***★★****★★**★*★★★**•★ 



INPUT PROCEDURE 

**•*•*•★**★★★**********★*★**********************************★★*•**★★★★*★**★★★★ * j 



void input(void)( 

clrscr( ); _setcursortype(_NORMALCURSOR) ; 



cpr i ntf ( "\r\n THE LONGITUDINAL FIN OF TRAPEZOIDAL PROF I LE\r\n" ) ; 

cpr i nt f ("\r\n ======================================================= ===========\r\n") ; 

cprintfC Fin Depth "); depth = i_L(L); 

cprintfC Heat Dissipation "); heat = i_H(H); 

cprintfC Taper Ratio = "); scanf ("%f ",&taper_ratio); 

cprintfC Emissivity = "); scanf("%f", Sepsilon); 
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cprintfC" Density ( kg / dnT3 ) = " ); scanf ("%f" # &density); 

cprintfC" Base Temperature "); Tbase = i_T(T); 

cprintfC 1 Thermal Conductivity "); k = i _K ( K ) ; 

cprintfC" Environmental Parameter (0 to .8) = "); scanf ("%f" # &V); 

cprintf C"\n======== ==================================================== ======\r\n"); 



if ( epsilon > 1.0 ]] epsilon <= OK 

cprintf ("\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 

getchC); exit(0);} 

if ( V < 0 [ J V > 0 . 8) { 

cprintf ("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 

getchC ); exit CO);} 

if ( taper_ratio > .99 ]] taper_ratio <= OK 

cprintf ("\r\n\r\nFATAL ERROR: Taper Ratio is not between 0 and 1.0"); 

getchC); exitCO);} 

_setcursortypeC_NOCURSOR); 
return; > 

/* *********************************************************************** 

GLOBAL FUNCTIONS 

A************************************************************************* ^ j 



double hrCvoidK 

float Tav = CTb + Ta) / 2; 

returnCsigma * epsilon * powCTav, 3) - CK2 / C2 * Tav)));> 
double mCvoidK 

returnC hrC) / Ck * s i nC f in_taper ) ) ) ; > 

double uCdouble x){ 

returnC 2 * sqrtCx * mC) ) );} 

double A_Cvoid){ 

returnC ua * C II Cua)*K0Cub) + I 0Cub)*Kl Cua) ) );> 

double B_Cvoid){ 

return C ub*ub*C lOCub)*KOCua)- I 0Cua)*K0Cub) ) / C 2*k*b_width*depth*mC ) ) );> 
double C_CvoidK 

return C 2*k*a_width*depth*mC )*ub*C 1 1 Cub)*K1 Cua) - 1 1 Cua)*K1 Cub) )/Cua) );> 
double D_Cvoid){ 

returnC ub*C 1 0Cua)*K1 Cub)+ 1 1 Cub)*K0Cua) ) );} 
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j * * * 



************************** 



************************* 



********* 



OPTIMIZE Analysis Program 
Longitudinal Fin of Triangular Profile 

**************************************************************************** * j 



U include 


<stdio.h> 




# include 


<stdl ib.h> 




# include 


<math . h> 




# include 


<conio.h> 




M include 


"uni ts.h" 




^include 


"bessel ,h" 




#def ine 


Yes 1 




#def ine 


No 0 




#def ine 


N 99 




#def i ne 


sigma 5. 66961 e- 8 




char 


f i le_name [20] ; 




double 


length, edge width, 


base_width 


double 


A, B, C, D, Ta, Tb, 


qb=0, ua, i 


float 


V, K2, epsilon, k, 


Tbase, dens 


int 


T, L, M, S, H, d, K, i, X, 2, ’ 


void 


input(), resul ts(), 


compute( ) ; 


FILE 


*data_in, *out; 





PN, heat; 



/* *************************** 



******************** 



***************** 



MAIN PROGRAM 



******************* 



*************** 



********************* */ 



int main(){ 

data_in = fopenC'units.ext", "r"); 

fscanf (data_in,"%d %d %d %d %d %d %d %d %d %d %f" # &T l &L # &M,&S l &H l &d,&K,&X l &Z l &Y f &error); 
fc lose(data_in); 
error = error / 100; 

input C ); 

heat = heat / (1 - V); 

Ap = (1.1234 * pow(heat/depth , 3)) / ( k * pow(si gma*epsi Ion, 2) * pow(Tbase,9) ); 

base_width = (2.4852 * pow(heat/depth # 2)) / ( k * sigma*epsi lon*pow( Tbase, 5 ) ) ; 

compute( ) ; 
resul ts( ); 
ex i t ( 0 ) ; 
return 0;} 

j* *********************************************************************** 



COMPUTE ROUTINE 

************************************************************************** * j 



void compute(K 



double 

double 

double 

double 

int 



A_( ) , B_(), C_(), D_( ) , m( ) # hr(), u(double x); 
length_inc, bl, db; 

inc, Tave, overall_diff = 100, Tmax=0, Tbv_min = 100000, data [N+1 ] [6] ; 
Al , B1, Cl, D1, Ttip, temp, qa, segmented i f f ; 
overflow, possible_overf l ow = No, number_of_tries = 0; 



length = 2 * Ap / base_width; 
edge_width = base_width / 100; 
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db = length / ( N + 1 ); 

bl = length; 

length = length - db; 

fin_taper = atan2( base_width, 2*bl ); 

K2 = 2 * epsilon * sigma * V * pow( Tbase, 4 ) ; 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2) 

else 

Ttip = Tave + 0.0001; 
inc = .5 * Ttip; i=l; 

while ( f abs(overal l_di f f ) > ( Tbase * error) j \ i <= N K 

overflow = No; i = 1; Ta - Ttip; Tb = Ta; 
segment_diff = 10; 

while ( fabs(segment_di f f ) > ( Tb * error) / N && overflow == No){ 

ub = u(db); 

temp = Ta / I0(ub); 

C = 2 * k * edge_width * depth * pow(m(),2) * ll(ub) / (ub * I0(ub)) 
qb = C * temp; 

segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_di f f ; 

if (Tb > Tbase +50 | | Tb < 0){ 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbv_min ) Tbvjnin = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;)} 

if (overflow == No){ 

data[i][0] = qb; data[i][l] = Tb; data [ i 3 [2] = 1; 
data [i ] [3] = 0; data [i] [4] = C; data [ i ] [5] = 1; 

Ta = Tb; qa = qb; ++i;} 



a = bl - length; 

while (i <= N && overflow == No){ 
segment_diff = 10; 

a_width = edge_width * a / ( bl - length); 
b_width = edge_width *(a+db)/(bl- length); 

while ( fabs(segment_di f f ) > ( Tb * error) / N && overflow == NoX 

ua = u(a) ; ub = u(a + db); 

A = A_( ) ; B = B_( ) ; C = C_(); D = D_(); 
temp =A^Ta+B^ qa; 
qb =C*Ta+D* qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_di f f ; 

if (Tb > Tbase + 50 | ] Tb < OH 

overflow = Yes; possibl e_overf low = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = f abs( Tbvjnin - inc); 
if (Ttip < Tave) Ttip = Tave;)} 

if (overflow == No){ 

data [ i ] [03 = qb; data [ i 3 Cl 3 = Tb; data [ i ] [2] = A; 
data [ i ] [3] = B; data [i ] [4] = C; data [ i ] [5] = D; 

Ta = Tb; qa = qb; a = a + db; ++i;}} 
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if (number_of_tr i es > 29 K 

clrscrO; cprintf ("Fatal Error: Unable to Compute"); 

getchO; exit(O);) 

if (overflow == NoK 

if ( Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overf low == YesK 
inc = inc / 2; 

if (over a l l_di f f >0) Ttip = Ttip * inc; 
else Ttip = Ttip + inc;) 

else 

Ttip = Ttip - overal l_di f f ; 
number_of_tr ies++; 
gotoxy(27,24); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f " , number_of__t r i es, w_T ( T , Tt i p) , w_T (T , Tb) ) ; 

if (Ttip < Tave) Ttip = Tave;}) 

for (i = N - 1; i >= 1; i--)( 

Al = A * data [ i ] [2] + B * data[i][4]; B1 = A * data [i 3 [33 + B * data[i)[53; 

Cl = C * data [i ] [2] + D * data [ i ] [4] ; D1 = C * data [i ] [33 + D * data[i][5]; 

A = AT; B = 81 ; C=C1; D=D1;> 

length = length + db; 

if (S == DC 

out = f open( f i le_name, "w"); 
for (i = 1; \ <= N; i++) 

fprintf(out, "%e # %e, %e\n", w_L(d, i * length / N), w_T (T , data [ i 3 [1 ] ) , w_H(H, 

data [ i 3 [03 )); 

fclose(out );) 



Ta = Ttip; 

mass = Ap * depth * density * 1000; 

PN = 2 * sigma * epsilon * pow( L ength , 2) * pow(Tbase,3) / ( k*base_width ) ; 
eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase,4) ); 
return; } 

/* ********************************************************************** 



OUTPUT PROCEDURE 

***★★**★***★★***★*★***********★**★★★★★★*********★★**★★★★**★*★★**★★★★★*★★*★ ★ j 



void resul ts( void)C 



c l rscr( ); 

cprintf ("\a OPTIMUM PARAMETERS"); 

cprintf ("\r\n THE LONGITUDINAL FIN OF TRIANGULAR PROFILE"); 

cprintf ("\r\n ==================r======r========r===============================\ r \n") ; 

cprintf(" Density = %4.2f kg / dm A 3\r\n" ,densi ty) ; 

cprintf (" Fin Depth = "); p_L(L, depth); 

cprintf(" Total Mass = "); p__M(M, mass); 

cprintf(" Fin Length = "); p_L(L, length); 

cprintfC Emissivity = %4.2f\r\n", epsilon); 

cprintf (" Profile Area = "); p_A(X # Ap); 

cprintfC Profile Number = %4 .3f \r\n",PN); 

cprintfC Base Heat Flow = "); p_H(H, qb); 

cprintfC Fin Efficiency = %4 .3f \r\n",eta) ; 

cprintfC Fin Base Width = "); p_L(d, base_width); 

cprintfC Base Temperature = "); p_T(T # Tb); 

cprintfC Edge Temperature = "); p_T(T,Ta); 

cprintfC Thermal Conductivity = "); p_K(K, k); 

cprintfC Environmental Parameter = %4.2f\r\n", V); 

cprintfC ================================================================= =\r\n" ) ; 
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cprintf ( " 
cpri ntf ( " 
cpri ntf ( " 
cprintf (" 
getch( ); 
return; > 



Transmission Parameter Matrix \r\n\r\n“); 
A = %4.4f B = ", A); p_2(2,B); 

C = "); p_Y(Y,C); 

D = %4.4f", D); 



jk kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



INPUT PROCEDURE 



kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk * j 



void input ( vo idX 

c l rscr( ); _s et cur sortype(_NORMAL CURSOR ) ; 



cprintf ("\r\n THE LONGITUDINAL FIN OF TRIANGULAR PROFILE"); 

cprintf ("\r\n ============================================================ ======\r\ n ") 

cprintfC" Fin Depth "); depth = i _L ( L ) ; 

cprintf (" Heat Dissipation "); heat = i_H(H); 

cprintfC Emissivity = "); scanf ("%f " , &epsilon); 

cprintfC Density ( kg / dnT3 ) = "); scanf("%f" # Sdensity); 

cprintfC Base Temperature "); Tbase = i_T(T); 

cprintfC Thermal Conductivity "); k = i_K(K); 

cprintfC Environmental Parameter (0 to .8) = "); scanf C%f", &V); 

cprintf ("\n== ============================================================== ==\r\n" ); 

if ( epsilon > 1.0 \\ epsilon <= OH 



cprintf C\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 

getch(); exit(O);) 

if C V < 0 ] | V > 0.8H 

cprintf ("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8“); 

getch(); exit(0);) 

_se t cursor type(_N0CURS0R) ; 
return; > 

j k kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk'hkkkkkkkkk'kkkkkkkkkkkkkkkkkkkkkkkkkkkk 

GLOBAL FUNCTIONS 



kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk*kkkk**kk***kkkkk*kkkkkkkkkkkk*kkkkkkkkkkkk k j 



double hr(void)( 

float Tav = (Tb + Ta) / 2; 

return( sigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)));) 
double m(void)( 

return( hr() / (k * si n( f i n_taper) ) ); ) 

double u(double x)( 

return( 2 * sqrt(x * m( ) ) );} 

double A_(void)( 

return( ua * ( I 1(ua)*K0(ub) + I0(ub)*K1 (ua) ) );) 

double B_(void){ 

return C ub*ub*C I0(ub) 1 * r K0(ua) - l0(ua)*K0(ub) ) / ( 2*k*b_width*depth*m( ) ) );> 
double C_(void)( 

return C 2*k*a_width*depth w m( )*ub*( I1(ub)*K1(ua) - 1 1 (ua)*Kl Cub) )/(ua) );) 
double D_(void)( 

return( ub*( l0(ua)*K1 (ub)+I 1 (ub)*K0(ua) ) );) 
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/* 



**************** 



************ 



********************************** 



UNIT CONVERSION Program Header File 
******************************************* 



***************************** 



^include <stdio.h> 
^include <conio.h> 
^include <string.h> 
^include <stdlib.h> 
^include "units. h" 



float i_T ( int T){ 
float V; 
switch (T ){ 

case 0: 
case 1 : 
case 2: 
case 3: 
return( V) ; ) 



cprintf ("(deg K) = "); scanf("%f", &V); break; 
cprint f ("(deg C) = "); scanf("%f" ( &V); V = V + 273.15; break; 
cprintf ("(deg F) = "); scanf("%f", &V); V = 5 * ( V + 459.4 ) / 9; 
cprintf ("(deg R) = "); scanf ('W, &V); V = 5 * V / 9;} 



void p_T(int T, float V)( 
switch (T )( 

case 0 
case 1 
case 2 
case 3 
return; } 



cprintf ("%4.2f deg K\r\n", V); break; 

V = V - 273.15; cprintf ("%4.2f deg C\r\n", V); break; 

V = 9 * ( V - 255.222222) / 5; cprintf ("%4 ,2f deg F\r\n", V); 

V = 9 * V / 5; cprintf ("%4.2f deg R\r\n", V);} 



break; 



float w_T(int T # float V){ 
switch (T)( 
case 0 
case 1 
case 2 
case 3 
return(V); ) 



break; 

V = V * 2 73 -15; break; 

V = 9 * ( V - 255.222222) / 5; break; 

V = 9 * V / 5;) 



void p__Z( int 1 , 
switch ( Z K 

case 0: 
case 1 : 
return; } 



float VX 



cpri ntf ("%4 .4f deg C / watt\r\n", V); break; 

V = 0.1760831 * V ; cpr fntf ("%4 . 34 deg F hr / BTUNrXn", V);} 



float w_Z(int Z, float V)( 
switch (Z){ 

case 0: break; 

case 1: V = 0.1760831 * V;} 

return(V); } 

void p_Y(int Y, float VX 
switch (Y){ 

case 0: cpr i ntf ( "%4 .4f watt / deg C", V); break; 

case 1: V = V / 0.1760831 ; cprintf ("%4.4f BTU / hr deg F" f V);} 

return;} 

float w_Y(int Y # float V){ 
switch ( Y )( 

case 0: 
case 1: 
return(V);) 



break; 

V = V / 0.1760831;} 



float i_L(int 
float V; 
switch (L)( 



LX 



case 


0: 


cprintf (" 


(m) 


= "); 


scanf ("%f " , 


SV); 


break; 


case 


1: 


cprintf (*' 


(cm) 


= "); 


scanf ("%f " , 


&V); 


V = 


V / 100; break; 


case 


2: 


cprintf (" 


(ft) 


= "); 


scanf("7.f" f 


£V); 


V = 


0.3048 V; break; 


case 


3: 


cprintf (" 


(in) 


= M ); 


scanf ( "%f " , 


&V); 


V = 


0.0254 * V; break; 


case 


4: 


cpr i ntf (" 


(mm) 


= ■'); 


scanf ( "%f " , 


&V); 


V = 


V / 1000;} 



break 
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return(V);) 



void p_L(int L, float V){ 
switch (L){ 

case 0: cpr ;*ntf ( "%4 . 2f m\r\n", V); break; 

case 1: V = 100 * V; cpr i ntf ( "%4 .2f cm\r\n", V); break; 

case 2: V = V / 0.3048; cprintf ("%4.2f ft\r\n" # V); break; 

case 3: V = V / 0.0254; cprintf ("%4.2f in\r\n", V); break; 

case 4: V = 1000 * V; cprintf ("%4.2f mm\r\n", V);} 

return;} 



float w_L(int L, 
switch (L){ 

case 0: 
case 1 : 
case 2: 
case 3: 
case 4: 
return( V);} 



float V) { 
break; 

V = 100 * V; break; 

V = V / 0.3048; break; 

V = V / 0.0254; break; 

V = 1000 * V;} 



float i_A(int AH 
float V; 
switch (AH 



case 


0: 


cpr i ntf ( 


" (rrT2) 


= "); 


scanf ("%f ", 


&V); break; 






case 


1: 


cpr i ntf ( 


"(cm" 2) 


= "); 


scanf ("%f", 


&V); V = V 


/ 


10000; break; 


case 


2: 


cpri ntf ( 


M ( f t *2) 


= "); 


scanf (“%f ", 


&V); V = V 


/ 


10.76; break; 


case 


3: 


cprintf ( 


M ( i n A 2) 


= "); 


scanf ("%f" # 


&V); V = V 


/ 


1550; break; 


case 


4: 


cpr i ntf ( 


"(fTYTl"2) 


= "); 


scanf ("%f ", 


&V); V = V 


/ 


1000000;} 



return( V) ; } 



void p_A(fnt A, float V){ 
switch (AH 

case 0: cprintf ("%4. 2 m * 2\r\n", V); break; 

case 1: V = 10000 * V; cprintf ("%4.2f cm A 2\r\n", V); break; 

case 2: V = V * 10.76; cprintf ("%4.2f ft A 2\r\n", V); break; 

case 3: V = V * 1550; cprintf ("%4.2f in A 2\r\n", V); break; 

case 4: V = 1000000 * V; cpr i ntf ( "%4 . 2f mm A 2\r\n*' # V);} 

return; } 



float w_A( int A, 
switch ( A ) { 

case 0: 
case 1 : 
case 2: 
case 3: 
case 4: 
return(V);} 



float V)( 
break; 

V = 10000 * V; break; 

V = V * 10.76; break; 

V = V * 1550; break; 

V = 1000000 * V;} 



float i_M(int M){ 
float V; 
switch ( M){ 

case 0: 
case 1: 
case 2: 
case 3: 
return(V);} 



cprintfC (kg) 
cprintf(" (gm) 
cprintf ("( Ibm) 
cprintfC (oz) 



"); scanf ( "%f " , 
"); scanfC%f", 
"); scanf ("%f " , 
"); scanf("%f", 



&V); break; 

&V); V = V / 1000; break; 

&V); V = V / 2.204622; break; 
&V) ; V = V / 35.27396; break;} 



void p_M( int M, 
switch (M ){ 

case 0: 
case 1 : 
case 2: 
case 3: 
return;} 



float VX 

cprintf ("%4.2f kg\r\n", V); break; 

V = 1000 * V; cprintf ("%4.2f grams\r\n", V); break; 

V = 2.204622 * V; cprintf ("%4.2f lbm\r\n" ( V); break; 

V = 35.27396 * V; cpri nt f ("%4 . 2f oz\r\n" ( V);} 



float w_M(int M, 
switch (MX 

case 0: 
case 1 : 



float V) { 



break; 

V = 1000 * V; break; 
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case 2: V = 2.204622 * V; break; 

case 3: V = 35.27396 * V;} 

return(V);} 



float i_H(int H)< 
float V; 
switch ( H ) { 

case 0: 
case 1 : 
case 2: 
case 3: 
return(V);) 



cprintf (“(watts) = 
cprintfC (kw) 
cprintf (“(BTU/sec) 
cprintfC (BTU/hr) 



"); scanf ( “%f “ , &V); break; 

= “); scanf ( M %f“ ( &V); V = 1000 * V; break; 

= ■•); scanf ( “%f “ ( &V); V = V * 1054.3502; break; 

= “); scanf ( "%f“ , &V); V = V / 3.414425;) 



void p_H( i nt H, 
switch ( H ) { 

case 0: 
case 1 : 
case 2: 
case 3: 
return;) 



float V){ 

cprintf ("%4.2f watts\r\n“, V); break; 

V = V / 1000; cprintf (“%4.2f kw\r\n", V); break; 

V = V / 1054.3502; cprintf ("%4.2f BTU / sec\r\n“ ( V);break; 

V = 3.414425 * V; cprintf ("%4.2f BTU / hr\r\n“ ( V);) 



float w_H(int 
switch ( H ) { 

case 0: 
case 1: 
case 2: 
case 3: 
return(V);) 



float V){ 
break; 

V = V / 1000; break; 

V = V / 1054.3502; break; 

V = 3.414425 * V;) 



float i_K( int K)< 
float V; 
switch (K){ 

case 0: cprintf ("(watts / m deg K) = "); scanf("%f", &V); break; 

case 1: cpri ntf (“(BTU / ft hr deg F) = "); scanf("%f“, &V); V = V / 1.731;) 

return( V); ) 



void p_K(int K, float V){ 
switch (KK 

case 0: cprintf (“%4.2f watts / m deg K\r\n“, V); break; 

case 1: V = 1.731 * V; cpri ntf ( “%4 . 2f BTU / ft hr deg F\r\n“ # V);) 

return; ) 
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Jk ************************************************************************** 



BESSEL FUNCTION Program Header File 16 



jk A************************************************************************** 

U include <stdio.h> 

U include <math.h> 

#include "bessel.h" 

double K0(float x)C 

double y, z; 

double P1=-0. 57721566, P2=0. 4227842, P3=0 . 23069756, P4=0. 0348859, P5=0. 00262698, 

P6=0. 0001075, P7=0. 0000074; 

double 01=1.25331414, Q2=-0 .07832358, 03=0.02189568, Q4=-0 . 01062446, Q5=0 .00587872, 

Q6=-0. 00251546; 
double 07=0.00053208; 

if ( x < 2 ){ 

y = ( x * x ) / 4; 

z = ( - log(x/2)*I 0(x) )+(P1+y*(P2+y*(P3+y*(P4+y*(P5+y*(P6+y*P7) ) ))));} 

elseC 

y = 2 / x; 

z = (exp( - x )/sqrt (x) )»(Q1+y*(Q2+y*(Q3+y*(Q4+y*(Q5+y*(Q6+y*Q7) ) ))));) 
return(z); ) 

double KKfloat x)( 

double y, z; 

double Pl=1 .0, P2=0. 15443144, P3=- .67278579, P4=-0. 18156897, P5=-0. 01919402, P6=-0. 001 10404, 

P7=-0. 00004686; 

double 01=1.25331414, 02=0.23498619, Q3=-0 . 03655620, 04=0.01504268, Q5=-0 .00780353, 

06=0.00325614; 

double Q7=-0.00068245; 

if ( x < 2 ){ 

y = ( x * x ) / 4; 

z = ( log(x/2)*I1(x)) + ( 1/x)*(P1 + y*(P2+y*(P3+y*(P4+y*(P5+y*(P6+y*P7)))) ));) 

else( 

y = 2 / x; 

z = (exp( -x)/sqrt(x))*(Ql+y*(Q2+y*(Q3+y*(Q4+y*(Q5+y*(Q6+y*Q7)) ))));} 
return(z);} 



double IOCfloat x)( 
double y, z, ax; 

double PI =1 , P2=3. 5156229, P3=3 .0899424 , P4=1 . 2067492, P5=0. 2659732, P6=0 .0360768, 

P7=0. 0045813; 

double 01=0.39894228, 02=0.01328592, 03=0.00225319, Q4=-0 . 001 57565 , 05=0.00916281, 

Q6=-0. 02057706; 

double 07=0.02635537, Q8=-0. 0164 7633, Q9=0 . 00392377; 

if (fabs(x) < 3.75 K 

y = ( x / 3.75 ) * ( x / 3.75 ); 
z = Pl+y*(P2+y*(P3+y* (P4+y*(P5 + y*(P6+y*P7) ) ) ) ) ;) 

else{ 

ax = fabs(x); 
y = 3.75 / ax; 

z = (exp(ax)/sqrt(ax))*(Q1+y*(Q2+y*(Q3+y*(Q4+y*(Q5+y*(Q6+y*(Q7+y*(Q8+y*Q9) )))))));} 
return( z);> 

double I1(float x)( 

double y, z, ax; 

double P 1 =0 . 5 , P2=0. 87890594, P3=0. 51498869, P4=0. 1 5084934 , P5=0 .02658733, P6=0. 00301 532, 

P7=0. 0003241 1 ; 

double 01=0.39894228, 02= -0 .03988024 , Q3=-0 . 00362018, 04=0.00163801, 05= -0.0103 1555, 
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□6=0.02282967; 

double Q7=-0 . 02895312 , Q8=0 . 01 787654 f Q9=-0 . 00420059; 

i f (f abs(x) <3.75 ){ 

y - ( x / 3.75 ) * ( x / 3.75 ); 

z = x* (Pl+y* (P2+y* (P3+y* (P4+y* (P5+y* (P6+y*P7) )))));} 

else{ 

ax = fabs(x); 
y = 3.75 / ax; 

z = (exp(ax )/sqrt (ax ))*(Ql+y*(Q2+y* (Q3+y* (Q4+y* (Q5+y* (Q6+y*(Q7+y*(Q8+y*Q9) ) ))))));> 
return(z) ;> 
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APPENDIX B COMPUTER CODE TO DERIVE OPTIMUM GEOMETRY 



!+ *************************************************************************** 
PROFILE AREA Analysis Program To Determine The Optimum Geometry 
Longitudinal Fin of Rectangular Profile 

**************************************************************************** ★ j 



# include <time.h> 

# include <stdio.h> 
#include <stdlib.h> 
# include <math.h> 
^include <conio.h> 



#def i ne 
#def i ne 
#def i ne 
#def ine 



Yes 1 

No 0 

N 100 

sigma 5.66961e-8 



double length, width_base, max_length; 

double A, B, C, D, Ta, Tb, qb=0, Heat, qb_max, db, Continue; 

float V, K2, epsilon, k, Tbase, Ap, taper_ratio, depth, error=0.01; 

int i, j, number_of_passes, failures; 

void computeO; 



j* *********************************************************************** 



MAIN PROGRAM 



************************************************************************** * j 

int main(K 
float dq; 

FILE *out ; 
random i ze( ) ; 
error = error / 100; 
c l rscr( ) ; 

for ( j = 1 ; j <102; j++){ 

number_of_passes = 0; 

epsilon = (float) ( random(98) ); 

epsilon = .01 + epsilon / 100; 

depth = (float) (20 + random(200 ) ) ; 

depth = depth / 100; 

k = 10 + random(200); 

Tbase = 100 + random(500); 

Heat = (float) (25 + random( 1000) ); 

Heat = Heat / depth; 

V = (float) ( random(80)); 

V = .01 + V / 100; 

V = 0; 

taper_ratio = (float) ( random(98) ); 
taper_ratio = 0.01 + taper_ratio / 100; 
taper_ratio = 1; 



Continue = Yes; 
qb max = - 1 ; 
qb”- 0; 
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dq = 0; 

Ap = 1 .635 * pow(Heat,3) / ( k * pow(sigma,2) * pow(epsi Ion, 2) * pow(Tbase,9)); 

length = .87 * Heat / ( sigma * epsilon * pow(Tbase,4)); 

db = length / 20016; 

length = length ■ db; 

width_base = Ap / length; 

if (depth < 1 0* w i dth_base) Continue = No; 
if (length < 10*width_base) Continue = No; 

while (Continue == YesK 

length = length + db; 
width_base = Ap / length; 
number_of_passes++; 
failures = 0; 

compute( ); 

while (qb == 0 && Continue == YesK 
fai lures++; 

length = length + db / 2; 
width_base = Ap / length; 
compute ( ); 
if (fai lures == 3)( 

Continue = No; 
number_of_passes = 1;}} 



dq = qb - dq; 

cprintf ("\r\nwidth_base = %4.3f length = %4.6f heat = %4.6f dq = %4.6f", 
wi dth_base*1000, length, qb, dq); 
dq = qb; 



if (qb_max < qb ){ 

max_length = length; 
qb_max = qb;) 

else 

Continue = No;) 

if (number_of_passes > 2K 
length = max_length; 
width_base = Ap / length; 
compute( ) ; 

cprintf ("\r\nWri ting Rec # %d d = %4.3f b = %4.6f qb = %4.3f Tb = %4.1f Ap = %4.1f", 
j, width base, length, qb, Tbase, Ap); 
out = f open( "OP*T_RECT . MAT" , "a" ) ; 

fprintf(out, "%d, %e, %e, %e, %e, %e, %e, %e, %e, %e, %e, %e\r\n", j, 

depth, length, Ap, k, Tbase, epsilon, V, qb, taper_ratio, width_base, width_base ) 
fc lose(out ); ) 

elseC 

cpr intf ("\r\nIgnoring This Set Of Parameters"); 

j=j - 1 ;>> 

return 0;) 

jk kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

COMPUTE ROUTINE 



★A************************************************************************ k J 



void compute ( vo id)( 

double Tav, inc, overall_diff = 100, Tmax = 0, Tbv_min = 100000; 

double A, B, C, D, m, hr, Yo, Ttip, temp, qa, Tave, segment_dif f ; 

int overflow, possible_overf low = No, number_of_tries=0; 

width_base = Ap / length; 
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K2 = 2 * epsilon * sigma * V * pow( Tbase, 4) ; 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if ( V==0 ) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = .0001 + Tave; 
inc = .5 * Ttip; i = 1 ; 

while ( f abs(overal l_di f f ) > ( Tbase * error) j j i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * width__base * epsilon * pow(Ta, 4) * depth; 

while (i <= N && overflow == NoK 

segment_diff = 10; 

while ( f abs(segment_di f f ) > ( Tb * error) / N && overflow == NoK 
Tav = (Tb + Ta) / 2; 

hr = sigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)); 
m = sqrt( ( 2.0 * hr) / ( k * width_base) ); 

Yo = sqrt( 2.0 * hr * k * width_base) * depth; 

A = cosh( m * length / N ); 

B = sinh( m * length / N ) / Yo; 

C = sinh( m * length / N ) * Yo; 

D = A; 

temp =A*Ta+B* qa; 
qb = C * Ta + D * qa; 
segment_diff = ( temp * Tb) / 2; 

Tb = temp + segment_di f f ; 

if (Tb > Tbase +50 | j Tb < 0)C 

overflow = Yes; possi bl e_overf low = Yes; 
if (Ttip < Tbv_min ) Tbvjnin = Ttip; 
inc = (Tbvjnin - Tmax) / 2; 

Ttip = fabs( Tbvjnin - inc); 
if (Ttip < Tave) Ttip = Tave;)) 

if (overflow == NoK 

Ta = Tb; qa = qb; ++i;)) 



if (number_of_tr ies > 29 )C 

Ta=0; Tb=0;qb=0; return;) 

i f (overflow == NoK 

if (Tmax < Ttip) 

if (Tb < Tbase ) Tmax = Ttip; 

overalljdiff = (Tb - Tbase) / 2; 

if ( possible_overf low == YesK 
inc = inc / 2; 

if ( overal l_di f f > 0 ) Tt i p = T t i p - inc; 
else Ttip = Ttip + inc;) 

else 

Ttip = Ttip - overal l_di ff ; 
number j>f _t r i es++ ; 
if (Ttip < Tave) Ttip = Tave;)) 

return; ) 
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j-k *************************************************************************** 

PROFILE AREA Analysis Program To Determine The Optimum Geometry 
Longitudinal Fin of Trapezoidal Profile 

******************************** *********************** ********************* * j 



U include <time.h> 
^include <stdio.h> 
^include <stdl ib.h> 
^include <math.h> 
tfinclude <conio.h> 
^include ‘'bessel.h" 



#def i ne 
#def i ne 
#def ine 
#def i ne 



Yes 1 

No 0 

N 100 

sigma 5.66961e-8 



double max_length # length, depth, edge_width, base_width, Ta, Tb, qb=0, inc, Tave; 

float dq, V, qb_max, db, epsilon, taper_ratio, k, error=0.01, Tbase, K2, Heat, Ap; 

int i, Continue, number_of_passes , j, Exit, failures; 

void computeO; 



J* *********************************************************************** 



MAIN PROGRAM 



★ A*********************************************************************'*** * j 



int main(K 

FILE *out; 
randomi ze( ) ; 
error = error / 100; 
clrscrC ); 

for ( j = 1 ; j < 102; j++X 

number_of_passes = 0; 

epsilon = (float) ( random(98) ); 

epsilon = .01 + epsilon / 100; 

depth = (float) (20 + random(200) ) ; 

depth = depth / 100; 

k = 10 + random(200); 

Tbase = 100 + random(SOO); 

Heat = (float) (25 + random( 1000) ) ; 

Heat = Heat / depth; 

V = (float) ( random(80)); 

V = .01 + V / 100; 

V = 0; 

taper_ratio = (float) ( random(98) ) ; 
taper_ratio = 0.01 + taper_ratio / 100; 

Continue = Yes; 

Exit = No; 
qb_max = - 1 ; 
dq - 0; 
qb = 0; 

Ap = (.63077 * taper_ratio + 1.31723) * pow(Heat,3) / ( k * pow(sigma,2) * pow(epsi l on, 2) * 
pow(Tbase,9) ) ; 

length = .82 * (-.0315113 * taper_ratio + 1.0425113) * Heat / ( sigma * epsilon * 
pow(Tbase,4) ) ; 
db = length / 20016; 

length = length - db; 

base_width = 2 * Ap / ( length * ( taper_ratio + 1)); 
edge_width = taper_ratio * base_width; 
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if (depth < 10*base_width) Continue = No; 
if (length < 10*base_width) Continue = No; 

while (Continue == Yes)( 

length = length + db; 

base_width = 2 * Ap / ( length * ( taper_ratio + 1)); 
edge_width = taper_ratio * base_width; 

f ai lures = 0; 
number_of_passes++ ; 

computet ) ; 

while (qb == 0 && Continue == YesH 
fai lures++; 

length = length + db / 2; 

base_width = 2 * Ap / ( length * ( taper_ratio + 1)); 

edge_width = taper_ratio * base_width; 

computet ); 

if (fai lures == 3)( 

Continue = No; 
number_of_passes = 1;}} 



dq = qb - dq; 

cpr i ntf ("\r\nwidth = %4 . 3 f length = %4.6f heat = %4.6f dq = %4.6f", base_width*1000 # 
length, qb, dq); 
dq = qb; 



if (qb_max < qb ){ 

max_length = length; 
qb_max = qb;) 

else 

Continue = No; ) 

if (number_of_passes > 2K 
length = max_length; 

base_width = 2 * Ap / ( length * ( taper_ratio + 1)); 
edge_width = taper_ratio * base_width; 
computet ) ; 

cprintf ("\r\nWri ting Rec # %d d = %4.3f b = %4.6f qb = %4.3f Tb = %4.1f Ap = %4.1f M , 
j, base_width, length, qb, Tbase, Ap); 
out = f open ( "OPT TRAP .MAT'*, “a*'); 

fprintf(out, "%d, %e, %e, %e, %e, %e, %e, %e, %e, %e, %e, %e\r\n" , j, 

depth, length, Ap, k, Tbase, epsilon, V, qb, taper^ratio, base_width, edge_width ); 
fclose(out);} 

else( 

cprintf ("\r\nIgnoring This Set Of Parameters"); 

return 0;) 

void computet void){ 

double m, hr, n, ua, ub. A, B, C, D, fin_taper, I1UA, KOUB, I0UB, KlUA, KOUA, I0UA, K1UB; 

double qa, bl , dH, a, a_width, b_width, Ttip, temp, segment_di f f , HUB; 

double b, inc, Tav, Tave, overall_diff = 100, Tmax = 0, Tbvjnin = 100000; 

int overflow, possible_overf low = No, number_of_tr i es = 0; 

bl = length / ( 1 - taper_ratio ); 
fin_taper = atan2( base_width, 2 * bl ); 
dH = length / N; 

<2=2* epsilon * sigma * V * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if ( v==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 
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Ttip = Tave + 0.0001; 
inc = .5 * Ttip; i = 1; 

while C f abs(overal l_di f f ) > ( Tbase * error) jj i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * edge_width * epsilon * pow(Ta, 4) * depth; 

a = bl - length; 

while (i <= N && overflow == No){ 
segmentjdiff = 10; 

a_width = edge_width * a / ( bl - length); 
b_width = edge_width * ( a + dH ) / ( bl - length); 

while (fabs( segmented! ff ) > ( Tb * error) / N && overflow == No){ 

Tav = C Tb + Ta) / 2; 

hr = sigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)); 
n = hr / (k * sin(f in_taper)); 
ua = 2 * sqrt(a * n ); 
ub = 2 * sqrt( (a + dH) * n ); 

1 1 UA = 1 1 (ua) ; 

K0UB = K0(ub); 

1 0UB = 1 0 ( ub ) ; 

K1UA = K1(ua); 

K0UA = K0(ua); 

10UA = IO(ua); 

K1UB = K1 (ub); 

HUB = 1 1 (ub) ; 

A = ua * ( 1 1UA*K0UB + 10UB*KlUA); 

B = ub*ub* ( 1 0UB*K0UA - I0UA*K0UB) / ( 2*k*b_wi dth*depth*n) ; 
C = 2*k*a_width*depth*n*ub*( I 1UB*K1UA - I 1UA*KlUB)/(ua); 

D = ub* ( I 0UA*K 1UB + !1UB*K0UA); 
temp =A*Ta+B*qa; 
qb = C * Ta + D * qa; 
segmentjdi ff = ( temp - Tb) / 2; 

Tb = temp + segment_di f f ; 

if (Tb > Tbase + 50 | j Tb < OK 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

i f (overf low == No){ 

Ta = Tb; qa = qb; a = a + dH; ++i;}} 



if (number_of_tries > 29 ){ 

T a=0; Tb=0;qb=0; return; } 

if (overflow == No){ 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overf low == Yes){ 
inc = inc / 2; 

if (overa l l_di f f > 0 ) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overa l l_di f f ; 
number_of_tri es++; 
if (Ttip < Tave) Ttip = Tave;}} 



return; } 
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j * *************************************************************************** 



PROFILE AREA Analysis Program Tb Determine The Optimum Geometry 
Longitudinal Fin of Triangular Profile 

**************************************************************************** ★ y 



U include <time.h> 
tfinclude <stdio.h> 
tfinclude <stdlib.h> 
tfinclude <math.h> 
#include <conio.h> 
tfinclude "bessel.h" 



tfdefine Yes 1 

#define No 0 

tfdefine N 99 

tfdefine sigma 5. 66961 e- 8 

double length, depth, edge_width, base_width, max_length; 

double Ta ( Tb, qb=0; 

float V, qb_max, db, epsilon, k, error=0.01, taper_ratio, Tbase, K2, Heat, Ap; 

int i, Exit, Continue, number_of_passes , j, failures; 

void computeC); 

I * ****★★★******★■*■****★***★****★**★***★★*★★*★★★★★★*★★★*★★★★★★****★★★★★*★★* 



MAIN PROGRAM 



*■*★★*****★***•★*****★★*★****★*★★********★★★******•****★★★★★★*★*★*•***★*★★*★*★ * j 



int main(X 
float dq; 

FILE *out; 
randomi ze( ) ; 
error = error / 100; 
c l rscr( ) ; 

for ( j = 1 ; j < 102; j++)L 

number_of_passes = 0; 

epsilon = (float) ( random(98) ); 

epsilon = .01 + epsilon / 100; 

depth = (float) (20 + random(200) ) ; 

depth = depth / 100; 

k = 10 + random(200); 

Tbase = 100 + random(500); 

Heat = (float) (25 + random( 1 000) ) ; 

Heat = Heat / depth; 

V = (float) ( random(80)); 

V = .01 + V / 1 00; 

V = 0; 

taper_ratio = (float) ( random(98) ) ; 
taper_ratio = 0.01 + taper_ratio / 100; 
taper_ratio = 0; 

Exit = No; 

Continue = Yes; 
qb_max = - 1 ; 
qb"= 0; 
dq = 0; 

Ap = 1.304 * pow(Heat,3) / ( k * pow(sigma,2) * pow(epsi Ion, 2) * pow(Tbase,9)); 

length = .95 * Heat / ( sigma * epsilon * pow( Tbase, 4 ) ) ; 

db = length / 20016; 

length = length - db; 

base_width = 2 * Ap / length; 
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if (depth < 10*base_width ) Continue = No; 
if (length < 10*base_width ) Continue = No; 

while (Continue == Yes)( 

length = length + db; 
base_width = 2 * Ap / length; 

number_of _passes++ ; 
failures = 0; 

compute ( ); 

while (qb == 0 && Continue == Yes)C 
f a i lures++; 

length = length + db / 2; 
base_width = 2 * Ap / length; 
compute( ); 
i f ( fai lures == 3 H 
Continue = No; 
number_of_passes = 1;» 



dq = qb * dq; 

cprintf ( "\r\ndel ta = %4.3f length = %4.6f heat = %4.6f dq = %4.6f", base_width*1000, 
length, qb, dq); 
dq = qb; 



i f ( qb_max < qb )( 

max_ length = length; 
qb_max = qb;) 

else 

Continue = No; ) 

if (number_of_passes > 2K 
length = max_length; 
base_width = 2 * Ap / length; 
compute( ) ; 

cprintf ("\r\nWrit ing Rec U %d d = %4.3f b = %4.6f qb = %4.3f Tb = %4.1f Ap = %4.1f", 
j, base_width, length, qb, Tbase, Ap); 
out = f open ("OPT TR I .MAT", "a"); 
edge_width = 0; 

fprintf(out, "%d, %e, %e, %e, %e, %e, %e, %e, %e, %e, %e, %e\r\n", j, 

depth, length, Ap, k, Tbase, epsilon, V, qb, taper_ratio, base_width, edge_width ); 
fclose(out);} 

elseC 

cprintf ("\r\nl gnoring This Set Of Parameters"); 

j=j-l;» 

return 0;} 

void compute(voidK 



double 
double 
doubl e 
double 
int 



m, hr, A, B, C, D, I1UA, KOUB, I0UB, KlUA, KOUA, I0UA, K1UB, HUB; 

fin_taper, bl , db, a, ua, ub, a_width, b_width; 

inc, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 100000; 

Tav, Ttip, temp, qa, segment_di f f ; 

overflow, possible_overf low = No, number_of_tri es = 0; 



edge_width = base_width / 100; 

db = length / ( N + 1 ); 

bl = length; 

length = length - db; 

fin_taper = atan2( base_width, 2*b1 ); 



K2 = 2 * epsilon * sigma * V * pow(Tbase,4); 
Tave = Tbase * pow(V, 0.25) + 0.0001; 



if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 
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else 

Ttip = Tave + 0.0001; 
inc = .5 * Ttip; i = 1 ; 

while (fabs(overal l_di f f ) > ( Tbase * error) |] i <= N K 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 
segment_diff = 10; 

while ( f abs(segment_di f f ) > ( Tb * error) / N && overflow == NoK 

Tav = (Tb + Ta) / 2; 

hr = sigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)); 
m = hr / (k * sin(f in_taper)); 
ub = 2 * sqrt( db * m ); 

I0UB = I0(ub); 

KlUB = Kl(ub); 

HUB = IKub); 

C = 2*k*edge_width* depth*pow(m, 2)*I 1UB/(ub * 10UB); 



temp = Ta / I0UB; 
qb = C * temp; 

segment_diff = ( temp - Tb) / 2; 

Tb = temp + segmented!' ff; 

if (Tb > Tbase +50 ] | Tb < 0)< 

overflow = Yes; possible_overf low = Yes; 
if (Ttip < Tbv_min ) Tbvjnin = Ttip; 
inc = (Tbvjnin - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == NoK 

Ta = Tb; qa = qb; ++ i ; } 

a = bl - length; 

while (i <= N && overflow == NoK 
segment_diff = 10; 

a_width = edge_width * a / ( bl - length); 
b_width = edge_width *(a+db)/(b1- length); 

while ( f abs(segment_di f f ) > ( Tb * error) / N && overflow == NoK 

Tav = (Tb + Ta) / 2; 

hr = sigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)); 
m = hr / (k * si n(f in_taper) ) ; 
ua = 2 * sqrt(a * m ); 
ub = 2 * sqrt( (a + db) * m ); 

1 1 UA = 1 1 (ua) ; 

K0UB = K0(ub); 

10UB = I0(ub); 

KlUA = K1 (ua); 

K0UA = K0(ua); 

I0UA = I0(ua); 

K1UB = K1 (ub) ; 

HUB = 1 1 (ub); 

A = ua * ( I 1UA*K0UB + I0UB*K1UA); 

B = ub*ub*( I OUB^KOUA - I0UA*K0UB) / ( 2*k*b_wi dth*depth*m) 
C = 2 1k k 1fc a_wi dth^depth 1k m 1fc ub 1fc ( 1 1UB*K1UA - 1 1UA*KlUB)/(ua) ; 

D = ub*( I 0UA*K1UB + I1UB*K0UA); 

temp =A*Ta+B lfc qa; 
qb = C * Ta + D * qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_di f f ; 
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if (Tb > Tbase + 50 ]J Tb < Q){ 

overflow = Yes; poss i bl e_overf l ow = Yes 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbvjnin - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;)) 

if (overflow == No){ 

Ta = Tb; qa = qb; a = a + db; ++i;>> 

if (number_of_tries > 29) ( 

Ta=0; Tb=0; qb=0; return;) 

if (overflow == NoK 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overal l jdiff = (Tb - Tbase) / 2; 

if ( possible_overf low == Yes){ 
inc = inc / 2; 

if (overal l_di f f > 0 ) Ttip = Ttip * inc 
else Ttip = Ttip + inc;) 

else 

Ttip = Ttip - overal l _di ff; 
number_of_tri es++; 
if (Ttip < Tave) Ttip = Tave;)) 

length = length + db; 
return; ) 
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APPENDIX C MATLAB EVALUATION M FILES 



J-k *************************************************************************** 



MATLAB EVALUATION M FILE To Determine The Optimum Geometry 
Longitudinal Fin of Rectangular Profile 

**************************************************************************** * j 

load rect; 

D = rect; 

N = 100; 

sigma = 5. 66961 e- 8; 



L 


= D ( : 


: , 2 ) ; 


b 


= D ( : 


: , 3 ) ; 


A 


= DC 


: f 4); 


k 


= DC 


:,5); 


To 


= DC 


,6); 


e 


= DC 


jy, 


V 


= DC 


:,8); 


qo 


= DC 


:.9); 


TR 


= DC 


: ,10) 


db 


= DC 


: # 11 ) 


de 


= DC 


: ,12) 


I 


= [1: 


: 1 01 ] 



C3 = A.* k.* ( ( (si gma*e) . A 2) .* To. A 9).* (L./qo). A 3; 
Cl = b.* sigma.* e.* (To. "4).* (L./qo); 

C2 = db.* k.* sigma.* e.* (To. A 5).* (L./qo). A 2; 

dispCRECT V = O') 
dispC ') 

disp( ' Constant One: Polynomial p is:') 

p = polyf i t(TR, Cl, 1) 

f i tl = polyval (p, TR); 

disp('Mean Square Error:') 

std_err = sqrt(sum( (Cl - f i 1 1 ) . A 2)/(N- 1 ) ) 

pause 

dispC ') 

disp( 'Constant Two: Polynomial p is: 1 ) 

p = polyf i t (TR, C2, 1 ) 

f i 1 2 = polyval (p, TR); 

disp('Mean Square Error:') 

std_err = sqrt(sum( (C2- f i t2) . A 2)/(N- 1 ) ) 

pause 

disp(' ') 

disp( 'Constant Three: Polynomial p is:') 

p = polyf i t( TR f C3, 1) 

fit3 = polyval (p, TR); 

disp('Mean Square Error:') 

std_err = sqrt (sum( (C3- f i t3) . A 2)/(N- 1 ) ) 
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Jir ************************************************* 

MATLAB EVALUATION M FILE To Determine The Opti 
Longitudinal Fin of Trapezoidal Profile 
**********"****************************************** 



load trap; 

D = trap; 

N = 100; 

sigma = 5.66961e-8; 



L 


= 


0 ( : , 2 ) ; 


b 


= 


D( : 


,3); 


A 


= 


D<: 




k 


= 


DC, 5); 


To 


= 


DC 


,6); 


e 


= 


DC 


, 7 ); 


V 


= 


DC 


,8); 


qo 


= 


D(: 


#9); 


TR 


= 


DC 


,10) 


db 


= 


DC 


,11) 


de 


= 


DC 


,12) 


I 


= 


[1:101] 



C3 = A.* k.* ( ( (si gma*e) . * 2 ) .* To. "9).* (L./qo). A 3; 
Cl = b.* sigma.* e.* (To. "4).* (L./qo); 

C2 = db.* k.* sigma.* e.* (ToC5).* (L./qo). A 2; 

di sp( 1 TRAP V = O') 
dispC ') 

di sp( 1 Constant One: Polynominal p is: 1 ) 

p = polyf i t(TR, Cl , 1 ) 

fitl = polyval(p,T R ); 

disp('Mean Square Error:') 

std_err = sqrt(sum( (Cl-f itl ). A 2)/(N-1 )) 

pause 

disp(' ') 

di sp( 'Constant Two: Polynominal p is:') 

p = polyf i t(TR, C2, 2) 

f i t2 = polyval (p, TR ); 

dispCMean Square Error:') 

std_err = sqrt ( sum( (C2- f i t2) . A 2)/( N- 1 ) ) 

pause 

disp(‘ ') 

disp( 'Constant Three: Polynominal p is:') 

p = polyf i t(TR , C3, 1 ) 

fit3 = polyval (p,TR); 

dispCMean Square Error:') 

std_err = sqrt ( sum( (C3- f i t3 ) . A 2)/(N- 1 ) ) 

pause 



************************** 
mum Geometry 



************************ * j 
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j-k **★★*★**★*****★★**★*★*★**•**★*★**★****★*★***★*★★***★***★★*★*********★**★★*★★ 



MATLAB EVALUATION M FILE To Determine The Opt i 
Longitudinal Fin of Triangular Profile 

-k-k-kk-k-k'k'k'k-k'k'k-k'k'k-k-k'k'k'k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-kk-k-k-k-k-k-k-k-k-k 



l oad t r i ; 

D = tri; 

N = 100; 

sigma = 5. 66961 e- 8; 



L 


= DC 


:,2) 


b 


= DC 


:,3) 


A 


= DC 


:,4) 


k 


= DC 


,5) 


To 


= DC 


:,6) 


e 


= DC 


:,7) 


V 


= DC 


.8) 


qo 


= DC 


:,9) 


TR 


= DC 


: , 10) 


db 


= DC 


-ID 


de 


= D(; 


,12) 


I 


= [1: 
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C3 = A.* k.* ( ( (si gma*e) . A 2 ) .* To. A 9).* (L./qo). A 3; 
C1 = b.* sigma.* e.* (To. A 4).* (L./qo); 

C2 = db.* k.* sigma.* e.* (To. A 5).* (L./qo). A 2; 

disp( ' TR I V = O' ) 
d i sp< 1 ') 

disp( 1 Constant One: Polynominal p is:') 

p = polyfit(TR, Cl, 1) 

fitl = polyval(p,TR); 

dispC'Mean Square Error:') 

std_err = sqrt (sum( (Cl - f i 1 1 ) . A 2)/(N- 1 ) ) 

pause 

dispC ') 

disp( 'Constant Two: Polynominal p is:') 

p = polyf i t (TR, C2, 1) 

f i t2 = polyval (p, TR); 

disp('Mean Square Error:') 

std_err = sqrt (sum( (C2- f i t2) . A 2)/(N- 1 ) ) 

pause 

disp(' ') 

disp( 'Constant Three: Polynominal p is:') 

p = polyfit(TR, C3, 1) 

fit3 = polyval (p, TR); 

disp('Mean Square Error:') 

std_err = sqrt (sum( (C3- f i t3). A 2)/(N-1 )) 



mum Geometry 
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